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Abstract 

Some ideas/concepts in relativistic heavy ion collisions are dis- 
cussed. To a large extent, the discussions are non-comprehensive 
and non-rigorous. It is intended for fresh graduate students of Homi 
Bhabha National Institute, Kolkata Centre, who are intending to pur- 
sue career in theoretical /experimental high energy nuclear physics. 
Comments and criticisms will be appreciated. 
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1 Introduction 

Surprisingly, our diverse universe consists of a handful of 'elementary' or 
'fundamental' particles. In FigJH I have listed the presently known elemen- 
tary particles. These elementary particles can be classified as (i) matter 
particles, the fermions and (ii) mediator particles, the bosons. The hand- 
ful of fundamental particles, can interact only in four definite manner, (i) 
strong interaction (ii) electromagnetic interaction (iii) weak interaction and 
(iv) gravitational interaction. In tableiH I have listed the mediators of the 
interactions, also shown the relative strength of the interactions. All these 
particles interact gravitationally. 

Study of strong interaction is generally called nuclear physics. Histori- 
cally, nuclear physics started with Rutherford's discovery of 'Nucleus' in his 
celebrated gold foil experiment (1909). The term 'Nucleus' was coined by 
Robert Brown, the botanist, in 1831, describing the cell structure (alterna- 
tively, by Michael Faraday in 1844), from the latin word 'Nux' which means 
'nut'. The result of gold foil experiment was so bizarre at that time that 
Rutherford commented like this, "It was almost as if you fire a 15 inch shell 
into a piece of tissue paper and it came back and hit you". The concept of 
Atomic Nucleus was completed with James Chadwick's discovery of 'Neu- 
tron' in 1932. Indeed, one can say that proper Nuclear Physics started in 
1932 after the discovery of neutron. 

For a long time 'Atomic Nucleus' supposed to be composed of protons (a 
term possibly coined by Rutherford for hydrogen nucleus) and neutrons and 
they are supposed to interact strongly. In the mean time there was much 
progress in the understanding of electromagnetic (EM) interaction. It was 
recognised that EM interaction arises due to exchange of photons between 
two charged particles. In analogy to EM interaction, in 1934 Hideki Yukawa 
put forward the hypothesis that strong interaction between nucleons origi- 
nate from exchange of mesons. At that time mesons were not known. He 
made this bold conjecture to obtain a theory analogous to electromagnetic 
interaction, where a photon mediates the force. He was only 27 years old 
then. In 1937 pions were discovered and in 1949 Yukawa was awarded the 
Noble prize in Physics. However, in later years, with the advent of parti- 
cle accelerators, experimentalists discovered hundreds of particles (mesons 
and baryons) many of which can be thought to be mediators of the strong 
interaction. People then tried to characterize those particles, study their 
internal symmetry [internal symmetry refers to the fact that one generally 
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Figure 1: Three generation of matter particles in the standard model. The 
mediator particles (Gauge bosons) are also shown. 

find a family of particles called multiplet, all with same or nearly same mass. 
Each multiplet can be looked upon as a realisation of some internal sym- 
metry]. I will not go into detail, suffice to say that Murray Gell-Mann and 
George Zweig (1964) found that all these particles, including protons and 
neutrons, consists of only a few building blocks which he termed as quarks. 
Murray Gellman picked the word 'quark' from the sentence 'Three quarks for 
Muster Mark' in James Joyce book, 'Finnegans Wake'. Simplest version of 
the quark model faces problem. Some baryons e.g. or then composes 
of identical quarks and violate Pauli's exclusion principle. To eliminate the 
contradiction, the concept of color was introduced. Color is a new quantum 
number. Only three colors required to be hypothesised. Murray Gell-Mann 
was born in September 1929. When he postulates quarks, he was 35 years 
old. He got Nobel prize in the year 1969. One can borrow G. H. Hardy's 
(known for discovering Ramanujan) words and say, 'creative physics is young 
man's game'. Take for example: Newton, at the age 23-24 gave the law of 
Gravitation, discovered Fluxions (calculus), Einstein discovered relativity at 
the age of 25-26. Wolfgang Pauli formulated his exclusion principle when he 
is 25 years old. 

Traditionally, nuclear physics is the study of nuclear matter at zero tem- 
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Table 1: Four fundamental forces, their relative strength and their mediators 
are listed. 
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perature and at densities of the order of the atomic nuclei, nucleon density, 
p ~ 0.17/m-^ or energy density e ~ 0.16 GeV / fm? . Advent of accelerators 
has extended the study to hundreds MeV of temperature and energy densi- 
ties several order of magnitude higher. At such high density/temperature, 
individual hadrons loss their identity and the matter is best described in 
terms of the constituents of the matter, e.g. quarks and gluons, commonly 
called Quark-Gluon-Plasma (QGP). Historically, T. D. Lee, in collaboration 
with G. C. Wick first speculated about an abnormal nuclear state, where 
mucleon mass is zero or near zero in an extended volume and non-zero out 
side the volume [l][2]- They also suggested that an effective way to search 
for these new objects is through high-energy heavy ion collisions. In this 
short lecture course, I will try to discuss some aspects of the matter at such 
high density and temperature. For a general introduction to the subject, see 
, [3] mis] [6] [7] i [9] . 



2 Conceptual basis for QGP formation 

For composite hadrons, with finite spatial extension, concept of hadronic 
matter appears to lose its meaning at sufficiently high density. Once we have 
a system of mutually interpenetrating hadrons, each quark will find in its 
vicinity, at a distance less than the hadron radius, a number of quarks. The 
situation is shown schematically in FigJ2J At low density, a particular quark 
in a hadron knows in partner quarks. However, at high density, when the 
hadrons starts to interpenetrate each other, a particular quark will not able 
to identify the quark which was its partner at lower density. Similar phe- 
nomena can happen at high temperature. As the temperature of a nuclear 
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Figure 2: Left panel shows a nucleus at normal density. The right panel 
shows the same at high density. 

matter is increased, more and more low mass hadrons (mostly pions) will 
be created. The system again will be dense enough and hadrons will starts 
to interpenetrate. The system where, hadrons interpenetrate is best consid- 
ered as a Quark matter, rather than made of hadrons. It is customary to 
call the quark matter as Quark-Gluon-Plasma (QGP). We define QGP as a 
thermalised, or near to thermalised state of quarks and gluons, where quarks 
and gluons are free to move over a nuclear volume rather than a nucleonic 
volume. Model calculations indicate that beyond a critical energy density 
~ 1 GeV/ fm^ , or temperature T^^ ~ 200 MeV, matter can exist only as 
QGP. 

QGP is the deconfined state of strongly interacting mater. Since at low 
density or low temperature quarks are confined within the hadrons and at 
high density or at high temperature, quarks are deconfined, one can talk 
about a confinement-deconfinement phase transition. I will discuss it later, 
but it turns out that the confinement-deconfinement transition is not a phase 
transition in thermodynamic sense (in thermodynamic phase transition, free 
energy or its derivative have singularity at the transition point), rather it 
is a smooth cross-over, from confinement to deconfinement or vice-versa. 
The mechanism of deconfinement is provided by the screening of the color 
charge. It is analogous to the Mott transition in atomic physics. In dense 
matter, the long range coulomb potential, which binds ions and electrons 
into electrically neutral atom, is partially screened due to presence of other 
charges, the potential become much more short range. 
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V{r) = el/r — )■ el/r x exp{—r/rci) (2.1) 

here r is the distance of the probe from the test charge Cq. td is the Debye 
screening radius and is inversely proportional to density, 

ro ~ n"^/^ (2.2) 

At sufficiently high density, r£, can be smaller than the atomic radius. A 
given electron can no longer feel the binding force of its ion, alternatively, 
at such density, coulomb potential can no longer bind electron and ion into 
a neutral atom. The insulating matter becomes a conducting matter. This 
is the Mott transition. We expect deconfinement to be the quantum chro- 
modynamic analog of Mott transition. Due to screening of color potential, 
quarks can not be bound into a hadron. Now one may wonder about the 
very different nature of QCD and QED forces. Interaction potential in QED 
and QCD can be expressed as. 



QED : V{r) ~ -e^/r (2.3) 
QCD : V{r) a/r + ar (2.4) 

While in QED, potential decreases continuously with increasing distance, 
in QCD, at large distance, it increases with distance. However, screening 
is a phenomenon at high density, or at short distance. The difference in 
QED and QCD at large distance is of no consequence then. More over, 
due to asymptotic freedom, in QCD interaction strength decreases at short 
distances, thereby enhancing the deconfinement. 

It may be noted that in insulating solid, at T > 0, conductivity is not 
exactly zero, it is exponentially small, 

aE ~ e-^^/^ (2.5) 

where AE is the ionisation potential. Above the Mott transition tempera- 
ture, o"^; is non-zero because Debye screening has globally dissolved coulomb 
binding between ion and electrons, but below the Mott transition tempera- 
ture, ionisation can produce locally free electrons, making aE small but non- 
zero. Corresponding phenomenon in QCD is the creation of quark-antiquark 
pairs in the form of a hadron. If we try to remove a quark from a hadron, the 
confining potential will rise with the distance of separation, until it reaches 
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Figure 3: Schematic representation of temporal evolution of universe in big 
bang theory. 



the value rriH, the lowest qq state. At this point, an additional hadron will 
form, whose anti-quark neutralises the quark we were trying to separate. 
This is the mechanism of quark fragmentation. 

2.1 Why QGP is important to study? 

QGP surely existed in very early universe. In FiglSl different stages of evo- 
lution of universe, in the Big bang model, are shown. 

(i) At the earliest time, temperatures are of the order of T ~ lO^^GeV, it 
is the Plank scale temperature. At this stage, quantum gravity is important. 
Despite an enormous effort by string theorists, little is understood about this 
era. 

(ii) We have better understanding of the later stage of evolution, say, 
around temperature T ~ 10^^ GeV. It is the Grand unification scale. Strong, 
and electroweak interactions are unified at this scale. The universe at this 
scale may also be supersymmetric (for each fermion a boson exists and vice- 
versa) . 

(iii) As the universe further expands and cools, strong and electroweak 
interactions are separated. At much lower temperature T ~ 100 GeV, elec- 
troweak symmetry breaking takes place. Baryon asymmetry may be pro- 
duced here. Universe exists as QGP, deconfined state of quarks and gluons 
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(iv) Somewhere around T ~ 100 MeV, deconfinement-confinement tran- 
sition occur, hadrons are formed. Relativistic Heavy Ion collider (RHIC) at 
Brookhaven National Laboratory (BNL), and Large Hadron Collider (LHC) 
at CERN, are designed to study matter around this temperature. 

(v) at temperature T ~ 1 MeV, nucleosynthesis starts and light elements 
are formed. This temperature range is well studied in nuclear physics ex- 
periments. For example at our centre (Variable Energy Cyclotron Centre, 
Kolkata), nuclear collisions produces matter around this temperature. 

(vi) at temperature T ~ 1 ev, universe changes from ionised gas to a gas 
of neutral atoms and structures begin to form. 

QGP may also exist at the core of a neutron star. Neutron stars are 
remnants of gravitational collapse of massive stars. They are small objects, 
radius ~ 10 Km, but very dense, central density ~ 10 normal nuclear matter 
density. At such high density hadrons loss their identity and matter is likely 
to be in the form of QGP. One important difference between QGP at the 
early universe and that in neutron stars is the temperature. While in early 
universe, QGP is at temperature T ~100 MeV, at the core of the neutron 
star it is cold QGP, T ~ MeV. Hot and dense matter with energy density 
exceeding 1 GeV/fim? may also occur in supernova explosions, collisions 
between neutron stars or between black holes. 

3 Kinematics of HI collisions 

Our knowledge of universe is gained through experiments. Horizon of human 
mind and of science is increased by solving puzzles posed by new and newer 
experiments. It is thus appropriate that we discuss kinematics of heavy ion 
collisions, which is very relevant for experimentalists. 
Throughout the note, I have used natural units, 

{i)h = c = kB = l, (ii) Metricifif^'' = diag{l, -1, -1, -1). 

When we calculate some observable, the missing h, c and ks must be put 
into the equation taking into account the appropriate dimension of the ob- 
servable. We also use the Einstein's summation convention, repeated indices 
are summed over (unless otherwise stated). Thus, 

d.j'^ = dtf + d,r + dyjy + d,r 
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Figure 4: A space-time diagram for the evolution of matter produced in 
relativistic heavy ion coUisions. 

3.1 Space-time picture 

FigJUdepicts the colhsion of two nuclei in (t,z) plane. Two Lorentz contracted 
nuclei approaching each other with velocity of light and collide at (t=0,z=0). 
In the collision process a fireball is created. The fireball expands in space- 
time going through various processes till the created particles freeze-out. In 
relativistic mechanics, neither At nor Ax are invariant distance. Invariant 
distance is Ar^ = At"^ — Ax"^. Appropriate coordinates in a relativistic 
collision is then proper time and space-time rapidity, 



proper time : r = yW—z^ (3-1) 
1 t + z 

space-time rapidity://^ = - In ^ (3-2) 

Region of space-time for which = — 2;^ > is called time like region, 
= — 2^ < is called space-like region, t = z line is called lightlike 
(only light or massless particles can travel along this line). Space-like region 
is inaccessible to a physical particle, it need to travel faster than light. For 
a massive particle, with speed v < c, only accessible region is the time-like 
region. Particle production then occurs only in the time like region. Space- 
time rapidity {rjs) is properly defined in the time like region only. 77^ is positive 
and negative infinity along the beam direction, t = ±z. rjs is not defined in 
space-like region . 
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3.2 Lorentz transformation 

In relativistic nucleus-nucleus collisions it is convenient to use kinematic vari- 
ables which take simple form under Lorentz transformation for the change of 
frame of reference. For completeness, we briefly discuss Lorentz transforma- 
tion. 

If x'^ is the coordinate in one frame of reference, then in any other frame 
of reference the coordinates x"^ must satisfy, 

g^udx'^dx'^ = g^i^dx^dx" (3.3) 

or equivalently, 

doc dx , > 

^-^^'=^- ^'-'^ 

The transformation has the special property that speed of light is same 
in the two frame of reference, a light wave travels at the speed \dx/dt\ = 1. 
The transformation — x"^ = A^x'^ + a^, being an arbitrary constant, 
satisfying Eq J3.4l i.e, 

g^,A';A>', = gp, (3.5) 

is called a Poincare transformation. Lorentz transformation is the special 
case of Poincare transformation when = 0. The matrix form a group 
called Lorentz group. 

A general Lorentz transformation consists of rotation and translation. 
Lorentz transformation without rotation is called Lorentz boost. As an ex- 
ample, consider the Lorenz boost along the z direction by velocity /3. The 
transformation can be written as, 



t'\^(l -h\ ft 
z' 1-/37 7 J U 



(3.6) 



where, 7 = 1/a/1 — /3^is the Lorentz factor. 



3.3 Mandelstam variables 

In FigJH a two body collision process is shown. Two particles of momenta pi 
and p2 and masses mi and m2 scatter to particles of momenta and p4 and 
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Figure 5: Pictorial diagram ofA + B^C + D process. 

masses 7713 and 7714. The Lorentz-invariant Mandelstam variables are defined 
as, 

= ml + ml + 2E^E2-2p^.p2 (3.7) 
^ = = (Pi - Ps)^ = (P2 - PaY 

= ml + ml-2ElE3 + 2p^.ps (3.8) 

u = = (Pi - ^4)^ = {P2 - Pzf 

= m? + m^-2EiE4 + 2pi.p4 (3.9) 

They satisfies the constrain, 

s + t + u = + m2 + + 7714 (3.10) 

3.4 Rapidity variable: 

In relativistic energy, rapidity variable, defined as, 

V = itafi^ (3.11) 

= llnl±^i4^ = tanh-i = tanh-i(/3i) (3.12) 

2 l-pjE \E) ^' ' ^ ' 

is more appropriate than the longitudinal velocity (/3l = p^/E). Rapidity has 
the advantage that they are additive under a longitudinal boost. A particle 
with rapidity y in a given inertial frame has rapidity y + dy in a. frame which 
moves relative to the first frame with rapidity dy in the —z direction. One 
can see this from the addition formula of relativistic velocity (3i and ^2- The 
resultant velocity. 
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is also the addition formula for hyperbolic tangents, 

, 1 / , N tanh(?/i) + tanh(|/2) .„^.^ 
+ = l + tanh(,Otanh(,,) ^'''^^ 

The underlying reason is that Lorentz boost can be thought of as a hy- 
perbolic rotation of the coordinates in Minkowski space. In terms of rapidity 
variable, velocity and Lorentz factor can be written as. 



/3 = tanh(y) 
7 = cosh(?/), 

and the transformation in Eq i3.6l can be rewritten as, 

/A /coshfe) -^Mft) (3.15) 

\z J y— smh[y) cosh(|/j J \z J 

which is a hyperbolic rotation. 

Rapidity is the relativistic analog of non-relativistic velocity. In the non- 
relativistic limit, p << m and Eq J3.11l can be written as. 



1 + + mVz 1 , ?Ti + mVz 
y = - In — = - In 



2 + ^2 _ 2 m-mvz 

= ]^[\Yi{l + Vz)-Hl-Vz)]^Vz (3.16) 

In terms of the rapidity variables, particle 4-momenta can be parame- 
terised as, 

= {E,p^,py,pz) = {rriT cosh y,p^,py,mT smb. y) (3.17) 
with transverse mass {rriT), 



TTiT = Jm^ +p'^= Jm^ +p1+pI (3.18) 
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3.5 Pseudo-rapidity Variable: 

For a particle emitted at an angle 6 with respect to the beam axis, rapidity 
variable is, 



2 E-p^ 

_ 1 A/m^ + + p cos 6* {"ilQ) 
2 a/ rn^ + — p cos 6 
At very high energy, p » m,the mass can be neglected, 

1 p + p cos 9 

y = - In 77 

2 p — p cos y 

= -lntan^/2 = r/ (3.20) 

?7 is called pseudorapidity. Only angle 9 determine the pseudorapidity. It 
is a convenient parameter for experimentalists when details of the particle, 
e.g. mass, momentum etc. are not known, but only the angle of emission is 
known (for example in emulsion experiments). 

3.6 Light cone momentum: 

For a particle with 4-momentum p{po,Pt,Pz), forward and backward light 
cone variables are defined as, 

p+ = Po+Pz (3.21) 
p_ = po-Pz (3.22) 

It is apparent that for a particle traveling along the beam axis, forward 
light cone momentum is higher than for a particle traveling opposite to the 
beam axis. An important property of the light cone is that in case of a 
boost, light cone momentum is multiplied by a constant factor. It can be 
seen as follow, write the momentum in terms of rapidity variable, p^ = 
{rriT cosh y,p^,py,mT sinh y) , 

p^ = mrey (3.23) 
p_ = nire-y (3.24) 
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3.7 Invariant distribution: 

E 

longitudinal direction is, 



Let us show is Lorentz invariant. The differential of Lorentz boost in 



dpi = ^{dpz - PdE) = -fidpz - 13?''^^'' 



E " 

= ^^{E - Pp.) = ^E* (3.25) 

where we have used, E^ = + Pr + p1 =^ EdE = Pzdpz- Then dpz/E 
is Lorentz invariant. Since pt is Lorentz invariant, d^p/E is also Lorentz 
invariant. 

The Lorentz invariant differential yield is, 

^d^N _ ^ d^N _ d^N 2 
d^p d^pj-dpz d^prdy 

where the relation dp./ E = dy is used. Some times experimental results are 
given in terms of pseudorapidity. The transformation from {y,PT) to {r],pT) 
is the following, 



dN L dN p^^^ 



d-qdpT y cosh y dydpx 

3.8 Luminosity: 

The luminosity is an important parameter in collider experiments. The re- 
action rate in a collider is given by, 

R = aL (3.28) 

where, a is the interaction cross section and L is the luminosity (in cm^^s^^), 
defined as, 

L = fn^ (3.29) 

where, 

/ ^revolution frequency , 

Ni,N2 = number of particles in each bunch. 
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n =number of bunches in one beam in the storage ring, 
A =cross-sectional area of the beams. 

3.9 Collision centrality 

Nucleus is an extended object. Accordingly, depending upon the impact 
parameter of the collision, several types of collision can be defined, e.g. cen- 
tral collision when two nuclei collide head on, peripheral collision when only 
glancing interaction occur between the two nuclei. System created in a cen- 
tral collision can be qualitatively as well as quantitatively different from the 
system created in a peripheral collision. Different aspects of reaction dy- 
namics can be understood if heavy ion collisions are studied as a function 
of impact parameter. Impact parameter of a collision can not be measured 
experimentally. However, one can have one to one correspondence between 
impact parameter of the collision and some experimental observable, e.g. 
particle multiplicity, transverse energy {Ex = EiSinOi) etc. For example, 
one can safely assume that multiplicity or transverse energy is a monotonic 
function of the impact parameter. High multiplicity or transverse energy 
events are from central collisions and low multiplicity or low transverse en- 
ergy events are from peripheral collisions. One can then group the collisions 
according to multiplicity or transverse energy. 

It can be done quantitatively. Define a minimum bias collision where 
all possible collisions are allowed. In Figj6] charged particles multiplicity 
{Nch) in a minimum bias collision is shown schematically. Minimum bias 
yield can be cut into successive intervals starting from maximum value of 
multiplicity. First 5% of the high N^^ events corresponds to top 5% or 0-5% 
collision centrality. Similarly, first 10% of the high N^^ corresponds to 0-10% 
centrality. The overlap region between 0-5% and 0-10% corresponds to 5-10% 
centrality and so on. Similarly, centrality class can be defined by measuring 
the transverse energy. 

Instead of impact parameter, one often defines centrality in terms of num- 
ber of participating nucleons (the nucleons that undergo at least one inelastic 
collision) or in terms of binary nucleon collision number. These measures 
have one to one relationship with impact parameter and can be calculated 
in a Glauber model. 
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Figure 6: Schematic representation of multiplicity distribution in minimum 
bias nucleus-nucleus collision. 



3.9.1 Optical Glauber model 

Glauber model views AA collisions in terms of the individual interactions of 
constituent nucleons. It is assumed that at sufficient high energy, nucleons 
carry enough momentum and are undefiected as the nuclei pass through 
each other. It is also assumed that the nucleons move independently in the 
nucleus and size is large compared to NN interaction range. The hypothesis of 
independent linear trajectories of nucleons made it possible to obtain simple 
analytical expression for nuclear cross section, number of binary collisions, 
participant nucleons etc. Details of Glauber modeling of heavy ion collisions 
can be found in ^U\. Below, salient features of the model are described. 

In FigJT] collisions of two heavy nuclei at impact parameter b is shown. 
Consider the two flux tubes, (i) located at a displacement s from the centre of 
target nucleus and (ii) located at a displacement s — b from the centre of the 
projectile nucleus. During the collision, these two flux tube overlap. Now, 
for most of the nuclei, density distribution can be conveniently parameterised 
by a three parameter Fermi function. 



(3.30) 
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Figure 7: Nucleus-nucleus collisions as viewed in a Glauber model. 

where po is the nucleon density, R the radius, a the skin thickness, w measure 
the deviation from a spherical shape. In tablej^l for selected nuclei, these 
parameters are listed. 



Table 2: Parameters R, a and of a three parameter Fermi distribution for 
selected nuclei are given. 



Nucleus 


R (fm) 


a(fm) 


w (fm) 


WQ 


2.608 


0.513 


-0.51 




4.2 


0.596 


0.0 




6.38 


0.535 


0.0 


208 


6.62 


0.594 


0.0 


238 


6.81 


0.6 


0.0 



p(r) in Eq J3.30[ normalised to unity, can be interpreted as the probability 
to find a given nucleon at a position r(= x, y, z). Then, 



dzpA{s,z), 



(3.31) 



is the probability that a given nucleon in the nucleus A (say projectile) is 
at a transverse distance s. Similarly, Tb(s — b) = J dzp{s — h, z) is the 
probability that a given nucleon in the target nucleus B is at a transverse 
distance s — b. Then Ta{s)Tb{s — b) is the joint probability that in an 
impact parameter b collision, two nucleons in target and projectile are in 
the overlap region. One then define a overlap function, at impact parameter 
b, 



TAB{h)= / dhTA{s)TB{s-h) 



(3.32) 
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Overlap function is in unit of inverse area. We can interpret it as the ef- 
fective area with which a specific nucleon in A interact with a given nucleon 
at B. If is the inelastic cross section, then probability of an inelastic in- 
teraction is cr7V7v^AB(b). Now there can be AB interactions between nucleus 
A and B. Probability that at an impact parameter b there is n interaction 
is, 

P{n, b) = (^"^^^ [aNNTABih)ni - aNNTAB{h)f''-'' (3.33) 

The first term is the number of combinations for finding n collisions out 
of AB collisions, the 2nd term is the probability for having n collisions and 
the 3rd term is the probability that AB — n collisions do not occur. 

The total probability of an interaction between A and B is 

, AB 

^ = ^ P(n, b) = 1 - [1 - cr;v7vT^B(b)]^^ (3.34) 

n=l 

Total inelastic cross-section is, 



/■oo 

a.nel = / 2Txhdh{l-[l-aMNTAB{h)]'^'') 
Jo 

2'jTbdb{l — exp{—aNNTAB{b)) (3.35) 



Total number of binary collisions is, 

N^oiiib) = J2 ^^(^' b) = ABTAB{h)aMN (3.36) 

The number of nucleons in projectile and target that interacts is called 
participant nucleons or the wounded nucleons. One obtains. 



Kart{h) = AjdhTA{s){l-[l-aNNTB{b-s)f) 

+ B j dhTB{b - s)(l - [1 - a^NTAis)]^) (3.37) 

Glauber model calculation of binary collision number or participant num- 
ber is energy dependent through the inelastic NN cross section (Ttvat. It 
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Figure 8: A optical Glauber model calculation for the impact parameter 
dependence of number of participant nucleons (Npart) and number of binary 
collisions (Ncou) in -\/i^^=200 GeV Au+Au collision. Inelastic nucleon- 
nucleon cross section is (TAr7v=42 mb. 



is common to take, aNN ~ 30 mb at Super Proton Synchrotron (SPS; 
yijv^ ^20 GeV), 40 mb at Relativistic Heavy Ion CoUider (RHIC; y/s^j^ ^200 
GeV) and 70 mb at Large Hadron Collider (LHC; v^^^ >1000 GeV). For 
demonstration purpose, in FiglHl I have shown a Glauber model calculation 
for Npart and Ncou as a function of impact parameter in a/s^^=200 GeV 
Au+Au collision. One understands that there is a one-to-one correspondence 
between impact parameter b and participant number or collision number. 

3.9.2 Monte-Carlo Glauber model 

In Monte-Carlo Glauber model, individual nucleons are stochastically dis- 
tributed event-by-event and collision properties are calculated averaging over 
many events. Optical Glauber model and Monte-Carlo Glauber model give 
very close results for average quantities like binary collision number or partic- 
ipant numbers. However, in the quantities where fluctuations are important, 
e.g. participant eccentricity, the results are different. Monte-Carlo Glauber 
model calculations proceed as follows: (i) nucleons in the colliding nuclei 
are distributed randomly following the probability distribution p(r), (ii) an 
impact parameter is selected randomly from a distribution dN/dh oc 6, (iii) 
assuming the nuclei are moving in the straight line, two nuclei are collided. 
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(iv) if the transverse separation between two colliding nucleons are less than 
the 'ball diameter' D = \/o^vn/^, they are tagged as interacted, and a reg- 
ister, keeping the coordinates of the colliding nucleons is updated. More 
details about the model can be found in p!0].|ll]. 

4 QGP and hadronic resonance gas in the 
ideal gas limit 

The basic difference between quarks and gluons inside a hadron and quarks 
and gluons in QGP as existed in early universe or in neutron star or as pro- 
duced in high energy nuclear collisions, is that as opposed to the former, 
the later can be treated as a macroscopic system. A macroscopic system 
is generally characterised by some state variables, e.g. number density (n), 
pressure (p), energy density (e), temperature (T) etc. Dynamics of the sys- 
tem is then obtained in terms of these state variables. In kinetic theory 
this programme is realised by means of a statistical description, in terms 
of 'one-particle distribution function' and its transport equation. Later, I 
will discuss some aspects of relativistic kinetic theory. Here, following [12], I 
present some simple calculations for the number density, energy density and 
pressure of a macroscopic system of particles of mass m, chemical potential 
/i, and at temperature T, when particles follow Maxwell-Boltzmann, Bose 
and Fermi-Dirac distributions. 



4.1 Maxwell-Boltzmann distribution 

Maxwell-Boltzmann distribution function is, 

/(rt=pi)3«p(^) (4.1) 

The distribution is of fundamental importance. Bose as well as Fermi- 
Dirac distributions can always be written as an infinite sum of Boltzmann 
distribution. 
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^^^^ (27r)3 e(^-A^)/^ ± 1 



oo 

_L_^(^)n.l,-n(.-.)/T (42) 



The ± corresponds to Fermi and Bose distribution respectively. 
For Boltzmann distribution, the number density, energy density and pres- 
sure can be obtained as, 



^ ^ d^pEexp ( ^ ~ ^ ) (4.4) 



(27r)3 7 ^ V ^ 
P = 7 — TTT / dy-^exp (4.5) 



(27r)3 y ^3 E V ^ 
Let us introduce the dimensionless variables, z and r 



m E a/IpI^ + JTi^ 



\p\ = TVt^ - z^, \p\d\p\ = T\dT, (4.7) 
\p\^d\p\ = T^tVt^ - z^dr, (4.8) 

In terms of r and z, the number density can be written as. 



27r 3 



POO 

/ dT{T^ - z'y/'re-^ (4.9) 

J z 



Closed form expression can be given for n in terms of the modified Bessel 
function of the second kind [13], 



Kn{z) has another representation which can be obtained from the Eq j4.1UI 
by partial integration. 
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2"-Hn- 1)! 1 



^ ' (2n-2)! z""' 



(4.11) 



Modified Bessel function lias a nice recurrence relation. If Kq and Ki 
are known, all the others can be easily obtained. For completeness, the 
recurrence relation is noted below. 



From Eq J4.11l one easily obtain. 



and the number density in Eq J4.9l can be written in a closed form. 



n 



27r2 



z^Ko(z\ 



T3 



2vr2 Vf) \t) 



Similarly, the energy density can be obtained as. 



(4.12) 



(4.13) 



(4.14) 



rp4 

2^' 

Now, from Eq JCTH 



J dT{r^-z^f'\^e~^ 



dr{T' - z^fl\{T^ - z'} + z')e 



2\-T 



zKi{z) = / drir^ - z^/^e-^ 

J z 

and the final expression for energy density is. 



(4.15) 
(4.16) 



rp4 



3^/^[3^2^2(^) + ^'i^l(^)] 



m 



m 



m 
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The expression for pressure is similarly obtained, 



— K2[ — )eT (4.18) 



27r2 VT/ VT/ 

The expressions for n, e and p are simplified in the massless limit, z = 
m/T — 0, when one can used the asymptotic relation for the modified Bessel 
function, 

lim z^'Kniz) = 2"-^(n - 1)! (4.19) 



JNumber density:n = ^— ^ V j^j VT"/ ^"^ 

m./T^o = -^eT (4.20) 



energy density :e 



m/T-s>0 



rp4 
1 

^2 



m 



gT 



(4.21) 



pressure:/? = i^y-f) ^2 (^^j 

_ M 

m/T^O — 

One do notice that for massless gas, p = relation is obtained. Above 
equations implicitly assumed that the degeneracy factor g = 1. If the particle 
has degeneracy g, the expressions for n, e and p has to be multiplied by the 
same. 



(4.22) 
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4.2 Bose distribution 

We write the Bose distributions as an infinite sum of Boltzmann distributions 

^^^^ " (27r)3 e(^-^)/^ - 1 

= J^sT.^+r'e--^'-"'^^ (4.23) 

To obtain close expressions, we will need Riemann zeta function. Rie- 
mann zeta function is a function of complex variable s = [x + iy) and ex- 
pressed as the infinite series, 

^(,) = yl = l + l + l + l + ... (4.24) 

n=l 

one can compute. 



C(o) 


= -1/2 


C(i) 


= oo 


C(2) 


= vrVe? 


C(3) 


= 1.202 


C(4) 


= 7rV90 



One also note an important relation, between Riemann zeta function and 
Dirichlet eta function, 

ri{s) = J2 = (1 - 2^-^)C(^), (4.25) 

n=l 

giving, 



r/(0) = 


(-l)C(O) = 1/2 


r/(2) = 


2^^ ^ 2 6 


r/(3) = 


^C(3) = ^1.202 


r/(4) = 


8^^ ' 890 
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Riemann zeta function or more precisely, Riemann hypothesis played and 
continue to play an important part in the development of mathematical the- 
ory. Riemann zeta function have trivial and non-trivial zeros. It has zeros at 
the negative even integers. Riemann hypothesis states that non-trivial zeros 
of zeta function has real part |, i.e. non-trivial zeros lie on the line ^ + it, t 
being a real number. The hypothesis is one of the most challenging problems 
in mathematics, and is not proved until now. Once Hilbert was asked about 
what would be in his mind if he is resurrected 1000 years later. He answered 
that he will inquire if Riemann hypothesis is proved. 

Let us now calculate the number density of a Bose gas. 



n 



dr'p 



exp \ ^ 1 - 1 



= p^/ rfV|:(+l)"'-p (!^) (4,26) 

If we define a temperature T' = T/n then above expression can be written 



as, 



n 



E 

n=l 



(27r) = 



d^pexp 



T 



(4.27) 



The bracketed quantity is just the Eq j4.3l which has been evaluated in 
Eq J4.14[ We immediately get. 



n 



1 /nrnx"^ /nm\ n 



n=l 



(4.28) 



Similarly, one obtain for energy density and pressure. 



p 



1 



n=l 



n=l 



(4.29) 



(4.30) 



27 



in the limit m — 0, — )■ 



E ;3 = ;;^C(3) = 1.202- (4.31) 



n 



n=l 



.3^ 3^ ,4.32) 



7r2 7r2 ' 30 

n=l 



E^A = ^C(4) = ^T' (4.33) 



TT^n^ 7r2'^ ' 90 

n=l 



4.3 Fermi distribution 

In analogy to Bose particles described in the previous section, for Fermion, 
number density, energy density, pressure can be written as. 



00 

n = 



n=l 

e = 



P = 

n=l 

In the limit m — 0, yU — )■ 



n=l 

(4.35) 

^ f:(-r'£;:.(^)^-^(^)- ("a) 

77,= 1 



n 



^E(-ir^'4 = ^lC(3) (4.37) 

n=l 

00 172 

2f^y(_l)"+i^ = l!LT4 (4.38) 

7^2 Z^v ; ^4 g3Q ^ ; 

n=l 

i^y(_l)-+iJ_ = l!LT4 (4.39) 

n=l 

One notes that in the massless limit, energy density, pressure in Bose and 
Fermi distribution differ by the factor | only. 



P 
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Table 3: Summary of properties of quarks 



quark 


symbol 


Charge 


constituent mass 


current mass 


flavor 




Q/e 


(MeV) 


(MeV) 


up 


u 


2/3 


~ 350 


1.7-3.1 


down 


d 


-1/3 


~ 350 


4.1-5.7 


strange 


s 


-1/3 


~ 550 


100^^° 


charm 


c 


2/3 


~ 1800 


1.29+[j-fi X 10^ 


bottom 


b 


-1/3 


~ 4 X 10^ 


4.191°:^^ X 10^ 


top 


t 


2/3 


~ 170 X 10^ 


172.9 ± 1.08 X 10^ 



4.4 Number density, energy density and pressure in 
QGP 

At high temperature, QCD coupling is weak and to a good approximation, 
quarks and gluons can be treated as interaction free particles. Gluons are 
massless boson, and Eqs J4.31|4.32l and 14.331 derived for massless bosons are 
applicable. However, they have to be multiplied by the degeneracy factor gg. 
For gluons, there are 8 colors and two helicity state and degeneracy factor is, 

dgiuon = color X spiu = 8x2 (4.40) 

Quarks are fermions with three color and two spin state. Also, for each 
quark, there is an anti-quark. Quarks comes in different (six in total) fla- 
vors. However, mass of all the quarks flavors are not the same. In tabl^ 
I have listed the constituent and current quark mass of the six known fla- 
vors. Current quark mass is the relevant mass here, it enters into the QCD 
Lagrangian. Constituent quark masses are used in modeling hadrons. In a 
sense they are dressed current quarks. As seen in tableJHl u and d quarks 
current mass is approximately same and can be assumed to be degenerate. If 
mass of Nf flavors are assumed to be same, the degeneracy factor for quarks 
can be obtained as, 

Qquark = pcLrticle — antiparticle x spin x color x flavor 

= 2x2 x3 X Nf (4.41) 

Considering that difference in distribution introduce an additional factor 
I in quark energy density/pressure, one can define a effective degeneracy 
factor for QGP, 
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7 

dQGP — dgluon + -^Qquark 
21 

= le + yiV/ (4.42) 

Now, in 1974, a group of physicist at MIT gave a model for hadron 
structure. The model become very popular and is known as MIT bag model 
[H]. In the model, the quarks are forced by a fixed external pressure to 
move only inside a fixed spatial region (bag). Inside the bag, they are quasi- 
free. Appropriate boundary conditions are imposed such that no quark can 
leave the bag. MIT bag model predict fairly accurate hadron masses. Color 
confinement is built in the model. However, chiral symmetry is explicitly 
broken at the bag surface. A remedy was suggested in cloudy bag model 

Equation of state (equation of state is a relation between the state vari- 
ables, pressure, energy density and number density) of QGP can be approx- 
imated by the Bag model. As in the bag model, in high temperature QGP, 
quarks are approximately free and even though it is a deconfined medium, it 
is confined in a limited region (albeit, confinement region is of nuclear size 
rather than of hadronic size). If B is the 'external bag pressure', the ex- 
pressions derived earlier can be augmented with the bag pressure to obtain 
energy density and pressure as. 



e = QQGp'^T' + B (4.43) 

P = ggop^T'-B (4.44) 

^ - 9QGP^T^ (4.45) 

In MIT bag model for hadrons, bag pressure B^^^ ~200 MeV. However, in 
the QGP equation of state, bag pressure is obtained by the consideration that 
QGP is a transient state and below a critical or (pseudo) critical temperature 
Tc, QGP transform into a hadronic matter or Hadron Resonance Gas. If the 
transformation is a first order phase transition, the Bag constant is obtained 
by demanding that at the transition temperature Tc, pressure of the two 
phases are equal. 



30 



PQGpiTc) = Phrg{Tc) (4.46) 

It will be discussed later, but explicit simulations of QCD on lattice in- 
dicate that for baryon free (/is = 0) matter, the transformation of QGP to 
HRG is not a phase transition in the thermodynamic sense, rather it is a 
smooth cross-over. In that case, thermodynamic variables in two phases can 
be joined smoothly to obtain the Bag pressure. 

4.5 Hadronic resonance gas 

QGP is a transient state. If formed in heavy ion collisions, it will cool back 
to hadronic matter at low temperature. At sufficiently low temperature, 
thermodynamics of a strongly interacting matter is dominated by pions. As 
the temperature increase, larger and larger fraction of available energy goes 
into excitation of more and more heavier resonances. For temperature T > 
150 MeV, heavy states dominate the energy density. However, densities of 
heavy particles are still small, pi ~ e~^^/'^ . There mutual interaction, being 
proportional to piPj ~ ^-{Mi+Mj)/T ^ suppressed. One can use Virial 
expansion to obtain an effective interaction. Virial expansion together with 
experimental phase shifts were used by Prakash and Venugopal to study 
thermodynamics of low temperature hadronic matter [16] . It was shown that 
interplay of attractive interactions (characterised by positive phase shifts) 
and repulsive interactions (characterised by negative phase shifts) is such 
that effectively, theory is interaction free. One can then consider interaction 
free resonances constitute the hadronic matter at low temperature. 

The expressions for energy density, pressure and number density for 
hadronic resonance gas, comprising hadrons, at temperature T and chem- 
ical potential p can be obtained by summing over the same for individual 
components of HRG, 



31 



AT 



£(T,/i) = ^£,(T,/i,) (4.47) 
1=1 

AT 

P(T,/.) = J]P.(2^,/^0 (4.48) 

i=l 

AT 

n(T,/i) = ^ni(T,/i,) (4.49) 



i=l 

The chemical potential fii is, 

/ij = Bifi + Sifis (4.50) 

where Bi = 0, ±1, ±2... and = 0, ±1, ±2... are the baryon and strangeness 
quantum number of the ith hadron. 

Earlier, I have derived the expressions for rii, Si andpj, for particles obey- 
ing Fermi distribution (Eqs J4.34|4.35[ 14.361) and for particles obeying Bose 
distribution (Eqs j4. 28114. 291 14.30p . They can be used in the above equations. 
However, in deriving those expressions it was implicitly assumed that par- 
ticles are point particles. The expressions can be corrected to account for 
finite size of hadrons. The correction is called 'excluded volume correction'. 
If Vi is the volume of the ith hadron, then available volume is, 

N 

V = V{l-Y,^,n,) (4.51) 

i=l 

One can estimate the excluded volume per particle as 1/2 of spherical 
volume of radius 2r/i, 

ISvr n 

Vi = v = —ri (4.52) 

Several procedures are in vogue to include the finite volume effect [IT] , [18] 
[19] , [20] , [21] , [22] . For example, in [17], [19] excluded volume effect is taken 
into account by reducing all the thermodynamic quantities including pressure 
by the reduction factor r = [1 + Y^. VjU^^iT, /^i)]~^- How ever the procedure 
is not thermodynamically consistent. Kapusta and Olive [20] advocated the 
following procedure, which is supposed to be 'thermodynamically' consistent. 
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Finite or excluded volume corrected pressure, energy density, temperature 
and entropy density are. 



Pxv - ^ ^ Ppt(T») (4-53) 

4t; 



1 



4ti 

1 , eptjT") 



(4.54) 



(4.55) 



where T* is the temperature of the system having point particles. B is the 
bag pressure, i?^/'*=340 MeV. 

In [18] the 'excluded volume model' pressure is expressed in terms of the 
ideal (point particle) gas pressure as, 

P(T, /i) = P^'^(T, j2),j2 = fi- vP{T, n) (4.56) 

For a given excluded volume f , Eq l4.56l can be solved to obtained pressure 
at a given temperature and chemical potential. Particle number density, 
energy density can be obtained as, 

e(T,/i) = ^""^^'^^ , (4.57) 
n{T,fi) = """"^^'^^^ , (4.58) 
See El], [22] for more details on excluded volume correction in HRG. 



5 Quantum chromo dynamics: theory of strong 
interaction 

Modern theory of strong interaction is Quantum Chromodynamics (QCD). 
Formally, QCD can be defined as a field theoretical scheme for describing 
strong interaction. QCD is built on three major concepts, (i) colored quarks, 
(ii) interaction between colored quarks results from exchange of spin 1 colored 
gluon fields and (iii) local gauge symmetry. 
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(i) Quarks: Quarks are fundamental constituents of matter. Quarks have 
various intrinsic properties, including electric charge, color charge, spin, and 
mass. Quarks can come in three colors (e.g. red, green and blue). In ta- 
bleJSl properties of the presently known quarks are listed. One notes that 
quarks posses fractional charges. Fractional charges are not observed in iso- 
lation. Millikan's oil drop type experiments give negative result for fractional 
charges. The experimental fact that quarks (fractional charges) are not ob- 
served in isolation, was accommodated in the theory by postulating 'color 
confinement'. Due to color confinement, quarks are never found in isola- 
tion. Quarks combine to form physically observable, 'color neutral', particles; 
mesons (pion, kaon etc.) and hadrons (protons, neutrons etc.) From tableJHl 
one can identify protons as composite of (uud) and neutrons as composite of 
(ddu). It may be mentioned here that the mechanism of color confinement is 
not properly understood as yet. QCD Lagrangian is highly singular at small 
momentum (large distance limit). Numerical simulation of QCD on lattice 
does indicate confinement. 

(ii) Gluons: Gluons are the mediators of the strong interaction. They 
are mass less bosons (spin 1). Indeed, role of photons in QED is played by 
gluons in QCD. But unlike photons, which are not self-interacting, gluons 
are. There are eight types of gluons. This can be understood if we note that 
quarks (anti-quarks) can carry three color charges. They can be combined 
in 9 different ways, 1 (singlet) colorless state and 8 (octet) colored states 
(3x3 = 1 + 8). Gluons can not occur in a singlet state (color singlet states 
can not interact with colored states). Hence there can only be 8 types of 
gluons. 

(iii) Gauge theory: QCD is a gauge theory, i.e. Lagrangian is invari- 
ant under a continuous group of local transformations. The Gauge group 
corresponding to QCD is SU(3). Below, I briefly discuss Gauge theory and 
SU(3) symmetry group. More detailed exposition can be found in text books, 

e.g.mm- 

5.1 Gauge theory in brief 

QCD is based on the principle of local gauge symmetry of color interac- 
tion. Here, I briefly describe the procedure to obtain local gauge symmetric 
Lagrangian. 

Consider a complex scalar field (f){x), with Lagrangian density. 
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£o(</'(x), d^(t>{x)) = d,ct>*d^<p - V{<P(t>*) (5.1) 
The Lagrangian is invariant under a constant phase change, 

^ f/ = e-^" (5.2) 

where a is an arbitrary real constant. This transformation is called 'global 
gauge transformation'. The theory is said to be invariant under global gauge 
transformation under the group U{\). Note U = e~*" is a unitary matrix in 
one dimension, UW = 1. The transformation, 

(j)' -> e-^"0, (5.3) 

is a global gauge transformation under U(l). 
If the complex field is written as, 

<P = ^(01+^02) (5.4) 

0* = ^(01-^02) (5.5) 
the transformation: — )■ e~*°</), (p* — )■ e*"0* gives. 



z0'2) = e-*"(0i + z02) (5.6) 



-^(0;-^0'2) = e^"(0i-#2)), (5.7) 



which is equivalent to. 



(p' = (01 cos a + 02 sin a) (5.8) 
0*' = (—01 sin a + 02 cos a) (5.9) 

The transformation 0(x) — )■ e~*°0(x) can be thought of as a rotation 
in some internal space by an angle a. Thus U(l) group is isomorphic to 
0(2), the group of rotation in two dimensions. (In group theory, two groups 
are called isomorphic when there is one to one correspondence between the 
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group elements. Isomorphic groups have the same properties and need not 
be distinguished). 

In a global gauge transformation, (f){x) must be rotated by the same 
angle a in all space-time points. This is contrary to the spirit of relativity, 
according to which signal speed is limited by the velocity of light. Then 
without violating causality, (p in all the spatial positions can not be rotated 
by the same angle at the same time. This inconsistency is corrected in local 
gauge transformation, where freedom is given to chose the phase locally, the 
phase angle a become space-time dependent, 

0(x) ^ U{x)(f){x) = e-^"(^V(a;) (5.10) 
Under such a transformation, 

dy{x) U{x)d''(f){x) + (j){x)d^U{x) (5.11) 

and Lagrangian is not invariant under the gauge transformation. The 'un- 
derlined' term must be compensated. This can be done by introducing a 
gauge field /l'^(x), which under the local gauge transformation transform as. 

A^(x) ^ A^(x) + -a^a(x), (5.12) 
e 

and replacing the partial derivative (9^) to covariant derivative (-D/^) defined 
as, 

D''(f){x) = [9^ + ieA^'{x)](f){x) (5.13) 

While, Lagrangian is now invariant under local gauge transformation, it 
is not the same Lagrangian as before. A gauge field yl'^(x) is now present as 
an external field. To obtain a closed system, we need to add a kinetic energy 
term, to be constructed from and its derivatives. The only term which is 
invariant under the gauge transformation is, 

= d^A" - d^A^" (5.14) 

Thus grangian density for a closed dynamical system, 

invariant under local U{\) gauge transformation. 



(5.15) 
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Lagrangian in Eq J5.15l is essentially for QED, which is a local gauge theory 
with U(l) group symmetry. Symmetry group for QCD on the other hand is 
SU(3). In contrast to U(l), which is an abelian group (i.e. group elements 
commute), SU(3) is non-abelian (group elements do not commute). Non- 
abelian nature of SU(3) group introduces additional complications. 

5.2 Brief introduction to SU(3) 

Special unitary group SU (N) is a Lie group isomorphic to that of all iV x TV 
special unitary matrices, 



detU=l (5.16) 
U^U = 1 (5.17) 

In general N x N complex matrices has 2N'^ arbitrary real parameters. 
The condition WU = 1 imposes N"^ condition and detU = 1 one condition. 
Hence SU{N) has A^^ — 1 arbitrary parameters. Correspondingly SU (N) has 
N"^ — 1 generators. La, obeying, 

7= AT 

[La, Lp] = i fap-yLry (5.18) 
7=1 

fa/s-y are the 'antisymmetric' structure constants (changes sign for inter- 
change of consecutive indices, fai3-y=— f 0^13= f -yap- One immediately notes 
that SU{3) is non-abelian, generators or the group element do not commute 
(in an abelian' group, structure constants are zero and generator and the 
group elements commute). For reference purpose, structure constants for 
SU(3) are noted in tableJH 

Naturally SU{3) has 8 generators. La = |Aq,, (a = 1, 2. ..8). are (3 x 3) 
Gell-Mann matrices, they act on the (color) basis states 




(5.19) 



For the sake of completeness, I have listed the 8 Gell-Mann matrices. 
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Ai 



A7 



/ 


1 


\ 


h 





° 









/ 1 


U 


U 


(: 


-1 


















—i 





























—i 


)■ 


i 









Table 4: Structure constant for SU(3) 



ijk 


123 


147 


156 


246 


257 


345 


367 


458 


678 


fijk 


1 


1 
2 


1 
2 


1 
2 


1 
2 


1 
2 


1 
2 


2 


2 



One does note that Gell-Mann matrices are generalisation of Pauli ma- 
trices; 



1 

1 



0-2 







0-3 







Mathematically, quark fields transforms as the fundamental representa- 
tion of color group SU(3). An infinitesimal element of the group is repre- 
sented by the transformation, 



X 

S 



Sx 



e 2 



(5.20) 
(5.21) 



where Ua{a = 1,2. .8) are arbitrary infinitesimal real numbers. 
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5.3 Lattice QCD 

Schematically, QCD Lagrangian has the form, 



£ = -^Fr^;, + Yl - '9y^> - ^^^] (5-22) 

flavors 

with, 

F;, = d,A: - d^Al + gnAX^l (5.23) 

here is the Gluon Gauge field of color a (a=l,2,...8),and m is the 'bare' 
quark mass, jahc is the structure constant of the Group and the quark 
spinors, 

( i^r\ 

^=\^A (5-24) 

Though the Lagrangian looks simple, it is not possible to solve it analyt- 
ically. Only in the high momentum regime, it can be solved perturbatively. 
Perturbative approach however fails in the low momentum regime. The rea- 
son being the running of the coupling constant. 

Running coupling constant refiect the change in underlying force law, as 
the energy/momentum scale, at which physical processes occur, varies. As 
an example, an electron in short distance scale can appear to be composed of 
electron, positron and photons. The coupling constant has to be renormalised 
to incorporate the change as the scale of physical processes varies. In QED, 
effective coupling constant at the scale q can be written as, 

a.„(9=) = ^ = ^ (5.25) 

QED coupling increase as the momentum scale is increased. In the other 
words, effective electric charge becomes much larger at small distances. 
In QCD coupling constant, at the momentum scale q is, 

aJq"^) = ^ = ,.2 "''^^ \.r T (5.26) 

where Nf is the number of flavors and A ^200 MeV is the QCD scale parame- 
ter. The coupling constant thus increases as the momentum scale q decreases. 
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Figure 9: Running strong coupling constant. 

Perturbative expansion in terms of coupling constant will not converge. In 
Figini I have shown the experimentally measured values of a^. Experimental 
measurements agree closely with QCD predictions. 

One possible way to obtain equation of motion is to simulate QCD on a 
lattice, i.e. solve the Lagrangian numerically. In lattice simulation, the space- 
time is discretized to reduce the infinite degrees of freedom of 'Field variables' 
to a finite and (numerically) tractable number. One immediately notices that 
due to finite dimension of the lattice, Lorentz invariance is broken. Gauge 
invariance however, is kept explicitly, by parallel transportation of the gauge 
fields between adjacent lattice sites. In the continuum limit, lattice spacing 
a — 0, Lorentz invariance can be restored. 

In the following, I briefly discuss some aspects of lattice QCD. Lattice 
QCD is intimately related to Feynman's path integral formulation of Quan- 
tum mechanics. Below I briefly sketch the ideas behind the path integral 
method and parallel transport. For more informative exposure on lattice 
QCD, see [21] , [25] , [26] , [27] . 
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5.3.1 Path integral method 

Richard Feynman is one of the most celebrated physicists of twentieth cen- 
tury. Apart from the path integral formulation of Quantum mechanics, he 
made pioneering contributions in Quantum electrodynamics, superfluidity 
and particle physics. He invented the diagrammatic approach of QED (the 
Feynman diagrams). In 1965, Feynman, along with Julian Schwinger and 
Sin-Itiro Tomonaga, was awarded Nobel prize for their contributions in QED. 




Figure 10: Particle trajectories or paths between (xj,to) and {xf,tj^) for 
discrete time steps. 

Consider the propagation of a particle from position Xi at time to to 
position Xf at time t^. For a given trajectory, {x,x), the action is, 

/•tN 

S= dtL{x,x), (5.27) 

J to 

L{x,x) is the Lagrangian. Path integral method states that the transition 
probability from {xi,to) to {xf,tN) can be expressed as the weighted sum of 
all the possible paths or trajectories, 

{xf{tN)\xi{to)) exp[zS{x{t),x{t))] (5.28) 

paths{P) 

Let us discretize the time interval into N steps, — to = NAt. In FigJTOl 
the discretized paths are shown . One understands that for small enough 
time steps, any continuous path can be adequately traced. Now, one can 
sum the trajectories at a particular time step, say t„. 
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/oo 
rfx(t„)e'^^"^*")'*(*"^) (5.29) 

The procedure can be repeated for each time steps. In the hmit N ^ oo, 



. N-l 



{xf{tN)\xi{to)) ^ M j JJ e^^("(*")'*(*"» (5.30) 

n=l 

M is some normahsation. 

It is easy to extend the formalism to fields. Consider a one dimensional 
field (f){x,t). Again, consider the transition amplitude for the field (f)i{x,ti) 

to (j)f{x,tf), 

(UxMM^^tf)) - 5^ e^^«^->'^^'^<^-» (5.31) 

(f>p 

As before, we discretize the time intervals in N steps. Additionally, we 
discretize the space coordinates into steps. Note that space is infinite 
dimension. Thus, discretization can only be an approximation of infinite 
space. 



{Ux,U)\^f{x,tf))^Af lim / n nd0(x™,tOe*''"^<"'"'*""'^'''*'^'"'"'*"^" 

m=l n=l 

(5.32) 

It is convenient to make a Wick's rotation, t = —ir so that the space is 
Euclidean. Then, 

iS = i J dtd^xC = — J drd^xCE = —Se (5.33) 

In terms of the Euclidean action {Se), the transition probability can be 
written as. 



,U)\<pAx,tf)) ~ M lim / n TT d0(x„,,t„)e-^^«<^(^'"'*")>'^^'^^(^'"'*")>) 

Nx,N-^oo J J- J- J- J- 

m=l n=l 

~ A/" y ©[^Je-^^ (5.34) 
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'D[</)] is the shorthand notation of the integration measure, 

N-1 

V[4>]=Y[l[d<P{Xm.,tn) (5.35) 

m=l n=l 

Now in statistical mechanics, central problem is to compute the partition 
function, defined as, 

Z = ^(0(x)|e-^^|0(x)), (5.36) 

the summation is over all the possible state \(j){x)). /3 = 1/T is the inverse 
temperature. It can be rewritten as. 



</- 

= ^(</.(x,t = -z/3)|</.(x,t = 0)), (5.37) 

The partition function in statistical mechanics then corresponds to the 
path integral formulation for the transition probability. 

Z =U f V[(P\e~^^ (5.38) 



This is an important realisation. All the tools of statistical mechanics 
can be applied to field theory problems. Expectation value of any observable 
can be obtained as, 

(O) = i 1 P[0]e-^-(9 (5.39) 
5.3.2 Parallel transport 

One of the problems in general relativity is the derivative of a vector (or 
more generally, a tensor) quantity. In fiat space-time, derivative of a vector 
can be computed easily, 

V-(.) = „n,i:!(^L±ii^I!M (6.40) 

h-^Q h 
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Figure 11: In a parallel transport, tangent angle of a vector is changed in 
curved space-time, but not in a flat-space-time. 

However, in a curved space-time, since the metric tensor g^" depend on 
space, additional terms arise. This can be understood from FigHT] In flat 
space, a vector at x, when transported to x + h, the tangent angle remains 
the same. But, in a curved space, the tangent angle is changed. In gen- 
eral relativity, this is accommodated by defining covariant (or semicolon) 
derivative, 

d;,^ = d.V + r^;„V^" (5.41) 
where T^^ is the Christoffel symbol and is defined as, 

~ 2^ V dx'' dx^ a?" 

T'^^V°' in Eq J5.41 1 accounts for the change in the vector's coordinate rep- 
resentation during the transport {T^^ = is flat space time). 
Covariant derivative, 

D^'(t){x) = [d^ + ieA^'{x)](f){x) (5.43) 

defined in Eq l5.13[ is analogous to parallel transport, ieA^{x) is the change 
in the field's representation during transport from x^ to x'^ + dx'^. Then, 



(5.42) 



(Pix^ + dx^) = (f){x^) +dx''Df,(f){x^) 

= (f){x^) + dx^'id^ + ieA^)(f){x^') 

= dx^'d^(j) + [1 + ieA^dx>')](f){x^) (5.44) 

The first term in Eq J5.44l is essentially a translational term. The 2nd 
term containing the A'^ describe the transport of gauge field between two 
close points x and x + dx. For infinitesimal distance, the 2nd term in Eq j5.44l 
can be written as, e^^^'^'^^^ (j){x^) . By repeated application of infinitesimal 
transport, the current (gauged) value of phase of the wave function 0, at 
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the 4-dimentional space-time point y is related to its value at some reference 
point X by the parallel transport, 

_ gie/J'dC''A,(c)^(^)^ (5.45) 

the integration in the exponent goes along some path C^y that connects x 
and y. For the non-abelian gauge group SU(3), a quark can alter its color 
under parallel transport. Then for SU(3) gauge fields, the exponential or the 
phase factor is a 3x3 unitary matrix. An extension of the above equation can 
be written as, 

0(2/) = Pe'^'^''^^"^M(?)^-0(x) (5.46) 

The symbol P means path ordering. To construct the matrix of parallel 
transport at finite distance, one has to subdivide the path C^y into small 
parts and form ordered product of parallel transport along these small parts: 

U = p^^9^dc^A'^^(0L^ ^ -Q^^ ^ igdC^Al{C)L^) (5.47) 

C 

U is the path dependent representation of an element of the gauge group 
G (presently SU(3)). 

5.4 Lattice formulation of QCD 

Lattice is a regular set of space-time points. A schematic representation of 
a lattice in two dimensions is given is FigJT^ For our purpose, we define, 

(i) site (node): the lattice points, characterised by the coordinate x, gen- 
erally in unit of the lattice spacing. 

(ii) Link: shortest distance connecting two sites, characterised by coordi- 
nates and direction, 

(iii) plaquette: elementary square bounded by 4-links, characterised by 
coordinates and two directions. 

In general, one also imposes periodic boundary conditions for bosons, 
Xn+i = Xi and anti-periodic boundary condition for fermions, Xn+i = —Xi 

Lattice QCD simulations are computer intensive. Total number of de- 
grees of freedom is very large on lattice. The fermions are defined on the 
nodes (site), '^/'^(i), where a(= 1,2,3) is the color index and a(= 1,2,3,4) is 
the Dirac index. They are complex, requiring 24 real variables per node. One 
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Figure 12: Schematic representation of a lattice, in two dimensions. 



associates the gauge fields with the links, U^^j, where i and j are the neigh- 
bouring points and a,b are color indices. U is a unitary 3x3 matrix, a total 
of 9 complex variables times four possible directions, i.e. 72 real variables 
per node for the link variables. In total in each node we have {2ANf+72) 
variables. For two fiavor QCD, even a small lattice will deal with 7864320 
real variables. Effectively, one has to compute a 7864320 fold integration. 

The relation between the matrices U and the gauge field A'^ is the fol- 
lowing, 

U^{x) = [expizA^dx^l A^ = J2 (5-48) 

a 

Ufj_{x) is the SU(3) matrix attached to the lattice link connecting the sites 
at X and x + dx, in the direction /x. Inverse of the matrix connects the sites 
in the opposite direction, 

f/_^(x + dx) = U-\x) = Ulix) (5.49) 
In lattice QCD, one evaluates the partition function, 

Z = J d[U]d[tlj]d[tlj']e~^^^^'^'^'^ (5.50) 

where the action Se{U,iP,4'') = J d'^xC^ipyip' , A) and d[ip] = Yln'^n repre- 
sents all the possible paths. 
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Gauge invariance is explicitly maintained In lattice QCD. As mentioned 
earlier, quark fields are placed on the nodes and gauge fields are associated 
with the links. One then parallel transports the gauge fields from lattice site 
n to n + 1, maintaining gauge invariance. Gauge invariant objects are made 
from gauge links between quark and anti-quark or products of gauge fields 
in a closed loop. In FigfT^ simplest close loop of gauge field is shown. It is 
called plaquette, product of 4 links connecting 4 adjacent nodes. 

• • • • 

: □ : 

X 




p 



Figure 13: A plaquette on the lattice. The unit vectors, fi, v are two generic 
directions. 

Pp.v{x) = U^{x)Ul{x)Ul{x + dx)U,{x + dx) (5.51) 
Let us consider each term separately ( we have omitted the "for ease). 



U^,{x) ^ exp{iaA^{x + 

^ exp{ia[A^{x) + a/29^A^(x)]) 

[/+(x) ^ exp{-ia[A^{x) + a/2duA^{x)]) 

U^{x + ae^) ^ exp[-iaA^{x + a{e^ + e^/2))] 

^ exp[-ia{A^{x) + ad^A^^x) + a/2d^Af,{x))] 

Uu{x + aCfj) ^ exp{ia[A^{x) + adfj,A^{x) + a/2dyAy{x)]) 



Product of the links then gives. 
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+ a'^[A^{x),A^{x)]) 

^ Lima-^oexp{ia^F^^{x)) 

^ 1 + za^F^, - ay2F^,F^'' + . 



(5.52) 



The term F^^ vanishes when summed over the indices fi and u and one 
obtain. 



a'F.^F'^'' = 2[1-P^,{x)] 



(5.53) 



Now the pure gauge action in the continuum, in terms of the scaled field 
A^' ^Ai", 

9 ' 



s 



9' 



(fx 



4 A"^ a 



(5.54) 



Comparing above two equations, pure gauge action on the lattice can be 
written as, 



(5.55) 



ll,l/,X 



5.5 Fermions on lattice 

Adding quarks to lattice action needs additional effort. Quark fields are 
defined on the nodes. Quarks are fermions and obey Pauli exclusion princi- 
ple. Thus they have to be included as anticommuting Grassmann numbers. 
Grassmann numbers are mathematical construction such that they are anti- 
commuting. A collection of Grassmann variables 6i are independent elements 
of an algebra which contains the real numbers that anticommute with each 
other but commute with ordinary numbers x 



OiOj = -9j9i (5.56) 
OiX = x9i (5.57) 
9^ = (5.58) 
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One also note that the operation of integration and differentiation are 
identical in Grassmann algebra, 



Ide = (5.59) 

ede = 1 (5.60) 

Grassmann numbers can always be represented as matrices. In general, 
a Grassmann algebra on n generators can be represented by 2" x 2" square 
matrices. 

In continuous Euclidean space-time, a fermion field ipi^) has the action. 



d^xi/j{x){i-f^'Df' - m)ij{x) (5.61) 
On the lattice is translate into, 

SpiU, i^) = J2 H^)MiU, X, y)^iy) (5.62) 
x,y 

where M is the Dirac matrix, essentially lattice rendering of the Dirac oper- 
ator, Ip + m. The functional integral for the partition function then become, 

[dU] [diP] [d^/;]e-^«(^)-^^(^''^) (5.63) 

Computing numerically with Grassmann variables is non-trivial. One 
generally integrate out the fermion fields, leaving only the gauge fields, 
weighted by the determinant of the Dirac matrix M, 

[dU]e-^^^^Uet[M{U)] (5.64) 

Before proceeding further, I must mention the well known problem of 
'Fermion doubling'. If fermion action is naively discretized on a lattice, 
spurious states appear. For each fermion on the lattice one obtain 2'^^^ = 16 
fermions. There are many ways to formulate Fermion action on a lattice, e.g. 
Wilson fermions, staggered fermions, domain wall fermions etc. We would 
not elaborate on them. We just mention that till today, Fermion action on 
lattice is inadequately treated. 



49 



5.6 Metropolis Algorithm 

Partition function in Eq J5.64l is a many fold integration. One generally uses 
Monte Carlo sampling to evaluate the partition function. One such algo- 
rithm is by Metropolis. There are several other algorithms also. Metropolis 
algorithm is based on the principle of detail balance. 
Metropolis algorithm proceeds as follows: 

(i) start from arbitrary configuration (e.g. randomly distributed), 

(ii) looks at the value of the field (say 0) at any given point and change 
it: 0', 

(iii) calculate the variation in action -.SS = S{(j)') — S{(f)). if 6S is negative, 
it is a lower energy state and desirable. One replace the old value with the 
new value 0'. If 6S is positive, one accepts the new value with the probability 
exp{—6S). 

The procedure, after many iterations will produce a equilibrium distribu- 
tion. Any physically relevant observable can be computed from the equilib- 
rium partition function, 

5.7 Wilson loop 

Consider a qq pair at a distance r. A schematic representation of the evolu- 
tion of the pair is shown in FigJTlk. In quantum mechanics, time evolution 
of the pair is governed by; 

^(t) = e-'^'i'ipit = 0) (5.66) 

For confining quark potential {V{r) ^ kr), as kinetic energy goes as 1/m, 
for infinitely heavy quarks. Egg ^ kr. In Euclidean space-time {t — —ir), 
time evolution of the pair is then governed by, 

where A is the area spanned by the qq system during its evolution. 

The Wilson loop is defined as the trace of the gauge fields along the world 
line. A typical Wilson loop is shown in FigJT^. It is just the product of link 
variables along the contour 
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(b) 



Figure 14: (a) Area spanned by the qq pair at a relative distance r, as a 
function of time, (b) a Wilson loop. 



w{r, T) = trUiU2...UN (5.68) 
In the continuum, expectation value of Wilson loop, for large T and r is. 



{w{r,T)) ~ (e 



-i dx'^Ai_. 



) 



-k.{Area) 



(5.69) 



The area law is a manifestation of confinement. 



5.8 Lattice QCD at finite temperature 

QCD at finite temperature can be simulated on a lattice where one of the 
4-dimension, say the time, is much smaller than the others. In the limit 
where the space dimensions go to infinity, but the time remains finite, the 
value of the temperature can be related to the time size. 



Time 



(5.70) 



Temperature 

Finite temperature QCD is then studied on a anisotropic lattice with, 



Nt « N, = Ny = 



(5.71) 



The central role in QCD at finite temperature is played the trace of the 
product Ux,^ along a line parallel to the time axis (see FigJTSl). The trace is 
called Polyakov loop. 

Consider two Polyakov loop separated by the distance d. Gauge invari- 
ance is ensured by periodicity of boundary condition which allows us to 
'close' the loops. The points denoted by A are physically same points due to 
boundary conditions. The correlation of the two loops as a function of their 
separation d decreases as. 




space 



Figure 15: Schematic representation of two Polyakov loop separated by a 
distance d. 



C{d) ~ e-^"^^'^^-' ~ e-^''('^)/^ (5.72) 

where Eqq{d) is the potential energy of the quark pair. Now imagine that 
one separates the two loop more and more such that one of the loop goes out 
of the lattice volume. Then one measure Eqq{d — oo), i.e. energy of a free 
quark. Therefore, the expectation value of 'one' Polyakov loop behaves as, 

(L) ~ e-^"?('^=°°)/^ (5.73) 

Polyakov loop can be identified as the order parameter of a confinement- 
deconfinement phase transition. 



confinement : Eqq{d = oo) = oo ^ (L) = 
deconfinement : Eqq{d = oo) = finite =^ (L) ^ 

(5.74) 

Now when ever there is a phase transition, some internal symmetry is 
broken. What is the symmetry broken in confinement-deconfinement phase 
transition? QCD has a hidden, discrete symmetry called 2'(3) symmetry. To 
understand the symmetry, let us define: 

Z{G): Centre of a group G is the set of elements that commute with 
every elements of G, 



Z{G) = \Z^G\zg = gz,g^a\ 



(5.75) 
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Figure 16: (left panel) Lattice simulations for energy density as a function 
of temperature, (right panel) Lattice simulations for pressure as a function 
of temperature. 



For SU(3), the center of group Z(3) has elements, (1, e*^'^/^, e*'^'^/^). One 
understand Z(3) symmetry as the group of discrete rotation around the 
unit circle in the complex plane. The Euclidean action is invariant under 
these group of rotation, but Polyakov loop is not. The issue of confinement- 
deconfinement is then related to breaking of Z(3) symmetry. In the confined 
phase {trL)=0 and Z{3) symmetry is preserved. In the deconfined phase, 
(trL) 7^0 and Z{3) symmetry is broken. 

5.9 Some results of lattice simulations for QCD equa- 
tion of state 

Several groups worldwide are involved in lattice simulations. Since these sim- 
ulations are costly, some groups have merged their resources to form bigger 
group. In the following I will discuss some representative results of lattice 
QCD. They are from Wuppertal-Budapest collaboration [28],|29]. However, 
similar results are obtained in simulation by other groups e.g. HotQCD |30] . 
As indicated above, in lattice QCD, one calculate the partition function, 

Z = f [dU]e-^^''^^^Y[det[M{U,mg)] (5.76) 

q 

where So is the Gauge action, f3 is related to the gauge coupling, (3 = l/g"^ 
and M is the Dirac matrix, rrig is the quark mass for flavor q. Once the par- 
tition function is known, all the thermodynamic variables can be calculated 
using the thermodynamic relations. 
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Figure 17: (left panel) Lattice simulations for entropy density as a function 
of temperature, (right panel) Lattice simulations for speed of sound as a 
function of temperature. 



In FigHni Wuppertal-Budapest simulations for energy density [e] and 
pressure (p), as a function of temperature is shown. One notes that e/T^ 
sharply rises over a narrow temperature range 150-200 MeV. At large tem- 
perature, it saturates. Very similar behavior is seen in simulated pressure, 
p/T^ saturates at large T. In FigJT6| the Stefan-Boltzmann limit is indi- 
cated. Simulated ejT^ as well as p/T^, though saturates, remains below the 
Stefan-Boltzmann limit. If we believe that at high temperature QCD matter 
exists as QGP, its constituents are not free, they are interacting. This is the 
reason QGP is call strongly interacting QGP (sQGP). 

In the left panel of FigjlTl Wuppertal-Budapest simulation for entropy 
density is shown. Entropy density over cube of the temperature also increases 
rapidly over a narrow temperature range T ^=150-200 MeV. At large tem- 
perature, it saturates below the Stefan-Boltzmann limit. e/T^, p/T^ or s/T^ 
are effectively proportional to the degeneracy of the medium. Temperature 
dependence of thermodynamic variable, e.g. energy density, pressure and en- 
tropy density thus indicate that effective degrees of freedom rapidly changes 
across the narrow temperature range T= 150-200 MeV. In the right panel of 
FigHT] variation of square of speed of sound (cs) with temperature is shown. 
Speed of sound shows a dip around temperature T ;^150 MeV. 

In FigJT8| renormalised Polyakov loop L^en on the lattice is shown. From 
a small value (~0) at low temperature, L^en increase at high temperature. 
However, the increase is not rapid, rather smooth and over a large inter- 
val of temperature. Smooth change of L^en indicates that the confinement- 
deconfinement phase transition is not a true phase transition, rather a cross- 
over. The cross-over temperature can be identified with the pseudo-critical 



54 




100 150 200 250 300 350 
T [MeV] 



Figure 18: Renormalised Polyakov loop on the lattice. 



temperature for the transition. It can be found by computing the inflection 
point of Lren ( SHI inflection point, curvature of a curve changes sign). For 
Wuppertal-Budapest simulation, cross-over temperature is ~ 160 MeV. 

5.10 Chiral phase transition 

We have talked about QCD conflnement-deconflnement phase transition. 
However, QCD has a well known phase transition called 'Chiral phase tran- 
sition'. Chirality means 'handedness'. Handedness can be understood from 
the helicity concept. Let us deflne the helicity operator, 

/i = J.p = (L + s)-p = s- p (5.77) 

h is the projection of spin on the momentum direction. For spin half fermions, 
helicity operator will have two eigen values, -1-1/2 and —1/2. A particle with 
helicity -1-1/2 (-1/2) is called right (left) handed particle. 

In FigJT9| two particles with helicity -fl/2 and -1/2 is shown. One un- 
derstands that for massive particles helicity is not a good quantum number. 
Massive particle will move with flnite speed v < c and one can go to frame 
from where particle will move backward and helicity will be reversed. How- 
ever, massless particles moves with speed c and helicity is a good quantum 
number for massless particles. 

Concept of chirality is more abstract. Consider a Dirac fleld ip for massless 
particle. The Lagrangian is. 
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Figure 19: Two particles with helicity +1/2 and -1/2 is shown. 

£ = i^Yd^ij (5.78) 
For the sake of completeness, we note that ip = V'Wo- We also list the 7 



matrices, 

.^(_°. o)-°^(o-0--("o)- 

and, 

75 = 75^ = ^7°7S'7' (5.80) 
7 matrices obey the anticommutation relations, 

{7^7n = 2^7'^^{7^7'} = 0, (5.81) 
Consider the following transformation, 

Av -.ip ^ e-'^^-^ = (5.82) 

a is the Pauli matrices and is the rotation angle. This is the general 
structure of a unitary transformation. The conjugate field transforms under 

Ay as. 

Ay : ^ ^ e+*?®^ = (1 + 2|e)^ (5.83) 



The Lagrangian is invariant under the transformation A 



iip^^p - iQ - (5.84) 

= (5.85) 

One say that the vector current V"^ = tjj'jf^^ip is conserved. 
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Let us now consider the following transformation, 



Ka-.^^ e-'^'^^ip = (l - ^75!©) ^ (5.86) 
=»V; ^ e-'^'^^ip = (l-i-f^^e^ij (5.87) 

where ant i- commutation relation 7075 = —7570 is used. The Lagrangian for 
massless Dirac particle transforms as, 



= (5.88) 

the 2nd term vanishes due to the anti-commutation relation {75, 7/^} = 0. 
The Lagrangian for massless Dirac particle is also invariant under the trans- 
formation Aa, with conserved 'Axial Current', = 'ijj'~ffj,'~f5^ip. 
Let us introduce the mass term in the free Dirac Lagrangian, 

6C = -rmpip (5.89) 
and see how it transforms under Ay and A^. 



Ay : mtpi) -> e"^*2®^e *2®^/; = m^'i' (5.90) 
A^ : miljilj = mipip — 2imQ{tlj—'~f^tlj) (5.91) 

Thus for massless Fermions, Dirac Lagrangian is invariant under the 
transformation, Aa_ and Ay, i.e. vector and axial vector currents are con- 
served. This symmetry is called Chiral symmetry and its group structure 
is SU{2)v X SU{2)a- For massive Dirac particles only the vector current is 
conserved. 

Chiral transition is signaled by the quark condensate (ipip)- In a chiral 
symmetric phase, (ipip) = 0. In the chiral symmetry broken phase (ipip) 7^ 0. 
In QCD, quarks masses are small but non-zero. Chiral symmetry is broken 
and quark condensate (tpip) 7^ 0. However, at sufficiently high temperature, 
quark mass decreases and condensate {ipip) ~^ ^iid one says that chiral 
symmetry is restored. 
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Figure 20: Subtracted Chiral condensate in lattice. 



Table 5: Some key properties of chiral and deconfinement transitions in 
QCD 





Chiral Phase 


deconfinement phase 




transition 


transition 


quark mass 





oo 


symmetry 


chiral symmetry 


Center group symmetry 


order parameter 


quark condensate 


Polyakov loop 



In Figj20l lattice simulation for quark condensate is shown. Generally to 
remove various uncertainties associated with lattice simulations, one calcu- 
late a subtracted quark condensate, 



s,T 



I = u,d 



(5.92) 



From the infiexion point of A; one computes chiral transition occur 
at Tc ^160 MeV. In Wuppertal-Budapest simulations, both confinement- 
deconfinement phase transition and chiral transition occur approximately at 
the same temperature However, the two transitions are unrelated. Some 
key properties of chiral transition and deconfinement transition is listed in 
tableJSl 
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Figure 21: Current understanding about nature of confinement- 
deconfinement phase transition as a function of quark mass m,j,mrf and m^. 

5.11 Nature of QCD phase transition 

In Figj2n the current understanding [31] about the nature of confinement- 
deconfinement phase transition, in a baryon free matter, as a function of 
quark mass m^j,mrf and is shown. The results can be summarised as 
follow: 

(i) In a pure gauge theory (m^ — oo), the transition is 1st order. 

(ii) For rriq — 0, the Lagrangian is chirally symmetric and there is a chiral 
symmetry restoration phase transition. It is also 1st order. 

(iii) For < rriq < oo, there is neither confinement- deconfinement phase 
transition nor a chiral symmetry restoring phase transition. The system 
undergoes a cross-over transition. The order parameter, e.g. Polyakov loop, 
or the susceptibility shows a sharp temperature dependence and it is possible 
to define a pseudo critical cross-over transition temperature. 

5.11.1 QCD phase diagram at finite baryon density 

At finite baryon density, Fermion determinant is complex and standard tech- 
nique of Monte-Carlo importance sampling fails. Several techniques have 
been suggested to circumvent the problem, (i) reweighting [32], [33] , (ii) ana- 
lytical continuation of imaginary chemical potential [51], [53] and (iii) Taylor 
expansion [3S],[SZ]- These methods has been used to locate the phase bound- 
ary in T — fiB plane, /i^ = 3/ig is the baryonic chemical potential. The 
calculations suggest that the curvature parameter in the expansion. 
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Figure 22: Lattice QCD calculation of QCD phase diagram, in imaginary 
chemical potential method [32] and Taylor expansion method [39] are shown. 
The red line is the chemical freeze-out curve obtained in a statistical model 



is small [38]. As an example, in Figj22| QCD phase diagram obtained in 
the analytical continuation method [32] (the filled circles) and in Taylor ex- 
pansion [39] (the filled squares) are shown. Both the methods gives nearly 
identical phase diagram for fiB/Tc{fiB = 0) < 3GeV, curvature parameter is 
small, K ~ 0.006. At larger fiB, they differ marginally. 

From theoretical considerations, QCD phase transition is expected to be 
1st order in baryon dense matter. Since at /ig ^ deconfinement transition is 
a cross-over, one expect a QCD critical end point (CEP) where the 1st order 
transition line ends up at the cross over. Location of the QCD critical end 
point is of current interest. At the critical end point, the first order transition 
becomes continuous, resulting in long range correlation and fluctuations at 
all length scales. Mathematically, it is true thermodynamic singularity. 

Experimental signature of QCD critical end point is tricky. Since at 
CEP, fluctuations exists at all length scale, one expects these fluctuations to 
percolate in the observables. Event-by-event fluctuations of baryon number, 
charge number can possibly signal a QCD CEP. 

In Figj22| chemical freeze-out curve [1Q],[1T], obtained in statistical model 
analysis of particle ratios are shown (the red line). Curvature of the chemical 
freeze-out curve is factor of 4 larger than the curvature in the QCD phase 



m 




(5.93) 
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diagram. Small curvature of the freeze-out curve, compared to the chemical 
freeze-out is interesting. Experimental signal of critical end point will get 
diluted as the deconfined medium produced at the critical end point will 
evolve longer to reach chemical freeze-out. Fluid will have more time to 
washout any signature of CEP. 



6 Color Glass Condensate 

In ultra-relativistic heavy ion collisions deconfined medium, called QGP can 
be produced. Theoretical considerations, however indicate that prior to 
QGP, a new form of matter, 'Color Glass Condensate (CGC) '[12], [13] may be 
formed. I briefly describe here the beautiful concept behind the color glass 
condensate. According to theory, the new form of matter (CGC) controls 
the high energy limit of the strong interaction [l2],[l3] and should describe, 
(i) high energy cross-sections, (ii) distribution of produced particles in high 
energy collisions, (iii) distribution of small x particles in a hadron and (iv) 
initial conditions for heavy ion collisions. 

As we know hadrons consist of gluons, quarks and anti-quarks. Con- 
stituents of hadrons e.g. quarks, gluons are generically called parton (the 
parton name was given by Feynman, while Murray Gellman picked the word 
'quark' from the sentence 'Three quarks for Muster Mark' in James Joyce 
book, 'Finnegans Wake'). At very high energy hadron wave function has con- 
tributions from partons, e.g. gluons, quarks and anti-quarks. A convenient 
variable to measure contribution of constituents to hadron wave function is 
the fraction of the momentum carried by the constituent (Bjorken x variable), 

-^constituent / Ej^adron (6'-^) 

Probability G{x) to obtain a parton with momentum fraction x and x + 
dx is generally called the parton distribution function. Parton distribution 
function depend weakly on the resolution scale Q^. One can write the density 
of small X partons as, 

^ = xG(x,Q2) (6.2) 
dy 

In Figj23| gluons distribution function as measured in HERA (Hadron 
Electron Ring Accelerator) is shown. One observes that gluon density rapidly 
increases at small x. It is also an increasing function of the resolution scale 



61 




10* 10-3 10-2 10-1 10" 



X 

Figure 23: gluon density measured in HERA in three momentum scale. 

(Q^). Increase in gluon density at small x is commonly referred to as the 
small X problem. It means that if we view the proton head on with increas- 
ing energy, gluon density grows. QCD is asymptotically free theory, coupling 
constant decreases at short distances. As the density increases, typical sep- 
aration between the gluons decreases, strong coupling constant gets weaker. 
Then higher the density, the gluons interact more weakly. However, den- 
sity can not be increased indefinitely, it will then lead to infinite scattering 
amplitude and violate the unitary bound (unitary bound is a constrain on 
quantum system, that sum of all possible outcome of evolution of a quantum 
system is unity). One then argues that as the gluon density increases repul- 
sive gluon interaction become important and in the balance, gluon density 
saturates. The saturation density will corresponds to a saturation momen- 
tum scale Qsat- Qualitatively, one can argue as follows: imagine a proton is 
being packed with fixed size gluons. Then after a certain saturation density 
or the closepack density, repulsive interaction will take over and no more 
gluon can be added to the proton. Naturally, the saturation density depend 
on the gluon size, for a smaller size gluon, the saturation density will increase. 
Then there is a characteristic momentum scale Qsat which corresponds to in- 
verse of the smallest size gluon which are close packed. Note that saturation 
scale only tells that gluon of size l/Qsat has stopped to grow. It does not 
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mean that number of gluons stopped to grow. 

It is very reasonable to assume that some effective potential V describe 
the system of gluons. If phase space density of gluons is denoted by p, 

_ 1 dN 

at low density, the system will wants to increase the density and V ~ —p. On 
the other hand , repulsive interaction balance the inclination to condensate, 
Vrepuision ~ ^sP^- These coutributious balance each other when p ~ 
Density scaling as inverse of interaction strength is characteristic of conden- 
sate phenomena such as super conductivity. 

Phase space density p = dyc^pr ~ ^/'^s can be integrated to obtain 
saturation momentum scale (Qs), 



-Ql (6.4) 



vri?^ dy as 

The origin of the name 'Color Glass condensate' is now clear. The word 
color refers to Gluons which are colored. The system is at very high density, 
hence the word condensate. The matter is of glassy nature. Glasses are 
disordered systems, which behave like liquid on long time scale and like solid 
on short time scale. The word 'glass' arise because the gluons evolve on time 
scale long compared to their natural time scale 1/Qsat- The small x gluons 
are produced from gluons at larger values of x. Their (the fast gluons) time 
scale is Lorentz diluted and can be approximated as a static fields. This 
scale is transferred to the small x gluons. The small x gluons then can be 
approximated as static classical fields. 

CGC acts as a infrared cut off when computing total multiplicity. For 
momentum scale pt > Qsati produced particles are incoherent and ordinary 
perturbation applies. For momentum scale pt < Qsat^ the produced particles 
are in a coherent state, which is color neutral on the length scale 1/Qsat- 

One may wonder about the quarks degrees of freedom. At high energy, 
gluon density grows faster than quark density and distribution is overwhelm- 
ingly gluonic. Fields associated with CGC can be treated as a classical fields. 
Since they arise from fast moving partons, they are plane polarised, with 
mutually orthogonal color magnetic and electric fields perpendicular to the 
direction of motion of the hadron. They are also random in two dimensions 
(see Fig|2i 
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Figure 24: CGC as high density gluon fields on a two dimensional sheet 
travelling near speed of light. 

There are many successful application of CGC model in explaining vari- 
ous experimental results. For completeness purpose, I will show two results 
obtained in CGC model [H]. In Figj25| in two panels, rapidity density of 
charged particles in pp collisions and energy dependence of charge multiplic- 
ity are shown. The solid lines in the figure are obtained in a CGC based 
model. It no small wonder, that CGC based model can explains the data. 
Such an description to the data, from a first principle model was not available 
earlier. 



7 Relativistic kinetic Theory 

QGP is a macroscopic system. Properties of many-body system depend on: 
(i) interaction of the constituent particles and (ii) external constraints. One 
characterises the system in terms of macroscopic state variables, e.g. parti- 
cle density, temperature etc. and of the characteristic microscopic parame- 
ters of the system. One then tries to understand certain equilibrium/non- 
equilibrium properties of the macroscopic system. In kinetic theory this 
programme is realised by means of a statistical description, in terms of 'one- 
particle distribution function' and its transport equation. From the transport 
equation, on the basis of conservation laws, hydrodynamic theory of perfect 
fiuid can be constructed. Supplementing the conservation laws with entropy 
law, hydrodynamics for dissipative fiuid is constructed. 

In the following, we briefiy discuss relativistic Boltzmann or the kinetic 
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Figure 25: (left panel) CGC model predictions for rapidity dependence 
dN/ drj of charged hadron multiplicities in proton - proton (antiproton) colli- 
sions as a function of the pseudorapidity at different energies. The data are 
taken from Ref . |45j . (right panel) Energy dependence of total multiplicity in 
proton - proton (antiproton) collisions. The vertical dotted line marks the 
LHC energies for proton-proton collisions (W = 14000GeV ). collisions. The 
experimental data are taken from Ref. |46j . 



equation. We then show that basic equations for hydrodynamics are obtained 
by coarse graining Boltzmann transport equations. Most of the discussions 
are from [T2] . 

7.1 Some basic definitions in kinetic theory 

(1) Distribution function, f{x,p): in kinetic theory, a macroscopic system is 
generally studied in terms of the distribution function, f{x,p). f{x,p)d^xd^p 
is defined as the average number of particles in small volume d^x, at time t, 
with momenta between p, p + dp. 

It is implicitly understood that particle content in the volume element 
d^x is large enough to apply concepts of statistical physics, yet, d^x small in 
macroscopic scale. 

(2) Particle four- flow A^'^: is defined as the 1st moment of the distribution 
function. 



65 



I ^P^'f{x,p) (7.1) 
4-components of particle 4-fiow can be identified as follows: 

Particle density: iV°(x) = J d^pf{x,p) (7.2) 

particle flow: A^*(x) = J d^p fix,p) 

= y"dV7(^,p),^ = 1,2,3 (7.3) 

where we have introduced the velocity u = p/p°. 

(3) Energy-momentum tensor T^*^: is the 2nd moment of the distribution 
function. 



T'%x) = j ^p^p-f{x,p) (7.4) 
The components can be identified as follows: 



energy density: T°°(x) = J d^pp^ f{x,p) 

energyflow: (.) = / #pA7(.,p), » = 1, 2. 3 

momentum flow or 

(4) Entropy four-flow S'^: 

S'-ix) = - j ^p^f{x,p)[\ogf{x,p) - 1] 



(7.5) 
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f{x,p) is a dimensionful quantity (dimension=/m~^Ge\^~^). To make 
it dimensionless, one generally multiply with and subtract unity. Note 
that absolute value of entropy is not measurable, only change in entropy is 
measurable. Then the observables remain unaffected. 

(5) Hydrodynamic four- velocity, u^: in each space-time point a time- like 
vector is defined, 

u>'{x)u^{x) = 1 (7.6) 

In the local rest frame, = (1, 0, 0, 0). 
with help of u'^ one defines a tensor quantity, 

A'^'ix) = g^"" -u^'u'' (7.7) 
It is called projector, annihilates that part of the 4-vector parallel to u^, 

A^^u^ = (7.8) 
Choices of hydrodynamic four-velocity: 

(a) Eckart's definition: Hydrodynamic four velocity is related to the par- 
ticle four flow A^'^, 

(7.9) 



(b) Landau and Lifshitz definition: u^^ is related to the flow of energy, 

u'^ = (7.10) 

In the study of high energy heavy ion collisions, central rapidity region is 
essentially particle free. It is difficult to define hydrodynamics four velocity 
according to Eckart's definition. Landau-Lifshitz choice of hydrodynamic 
velocity is preferred as it is related to energy flow. 



7.2 Physical quantities of a simple system 

With the help of hydrodynamic four velocity u{x) one can define relevant 
macroscopic quantities, e.g. particle density, energy density, heat flow, the 
pressure tensor and entropy density, in a covariant manner. 

(i) particle density is the density of particles in the rest frame m(1, 0, 0, 0), 
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n = N''u^ (7.11) 

(ii) energy density of the particles in the rest frame, 

e = u^T'^'u, (7.12) 

(iii) heat flow: the difference of energy flow and flow of enthalpy h carried 
by the particles, 

= (w^T'^- _ /iA^-)A/^ (7.13) 
where enthalpy per particle is defined as, 

h = ^ (7.14) 

n 

with p the local hydrostatic pressure, to be defined shortly. Heat flow has 
the property that it is transverse to the hydrodynamic 4-velocity, 

g% = (7.15) 

(iv) Pressure tensor: 



A^T^'A^ (7.16) 



It is symmetric when T'^" is symmetric. In the local rest frame, it is 
purely spatial, 

Pll = 0; p'lr = pL'r = 0;PI'r = ^, j = 1,2,3 (7.17) 

Pressure tensor has a 'reversible' and an 'irreversible' part, 

p is the hydrostatic pressure. The quantity 11'^'^ is called the viscous pressure 
tensor. Writing Eqs J7.12|7.13l and l7.18l in a slightly different manner, we find, 



e = u^T'^'u, (7.19) 
q^^ + hA^"'N^ = u^T'^'^A^ (7.20) 
-pA^"' + U^"' = A^r""A^ (7.21) 
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The terms can be collected to obtain the expression for energy-momentum 
tensor, 

T'"' = TZ + Ttrr, (7-22) 
with the 'reversible' and 'irreversible' parts, 

TZ = eu^u-" - pA"-" (7.23) 
TZ = [(g^ + /lA'^'^iV,)^" + (g'^ + /iA'^'^iV,)^'^] + n^'^ (7.24) 

As mentioned earlier, two choices of hydrodynamic velocity is popular. 
In the Eckart frame, = N^'/^/WN^, and A^'^N^ = 0. Heat flow is, 

q" = u^T'^'A^ (7.25) 
and the irreversible part of the energy-momentum tensor reduces to, 

Ttrr = [q^u" + q^u^] + W (7.26) 



In the Landau-Lifshitz frame, = T'^^'^u'^ / u^Tp'^T^t-u'^ . Using the 
property, A^'^u^, = 0, heat flow is, 

q^ = -hA^'^Nu (7.27) 

The irreversible part of energy-momentum tensor then has the simplified 
form. 



=W (7.28) 

In Landau-Lifshitz frame, the energy momentum tensor does not contain 
heat flow q^. However, as manifest in Eq J7.27[ heat flow is not zero. It is 
manifested in particle flow. Using the projector, A'^'^, the particle four flow 
can be split into two part, in the direction of four velocity and in the direction 
perpendicular to it. 



In the Landau-Lifshitz frame. 



(7.29) 
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h 



(7.30) 



It is interesting to note that in the Eckart's frame, oc A^'^, particle flow 
is in the direction u'^. While energy-momentum tensor explicitly contain heat 
flow, is identically zero. 

The irreversible part of the energy-momentum tensor Tj^'^ leads to dissi- 
pation. In a realistic system, Tj^'^ is never identically zero. However, if it is 
neglected, the fluid is called ideal or inviscid fluid. 



7.3 Relativistic Kinetic (TRANSPORT) Equation: 



In kinetic theory, macroscopic system is described in terms of one body distri- 
bution function f{x,p,t). Kinetic or transport equation give the space-time 
development of the distribution function. It was originally derived by Boltz- 
mann and called Boltzmann equation. Ludwig Eduard Boltzmann (1844- 
1906) was an Austrian Physicist. He made major contributions in statistical 
thermodynamics. He was an early proponent of atomic theory, when atomic 
theory has not gained ground. During his life time, his works were not ap- 
preciated. Famous physicists like Ernest Mach, Wilhelm Ostwald opposed 
his atomic view. Rejection of his views made Boltzmann depressed and in 
1906, he committed suicide. In the following, we derive the relativistic ver- 
sion of the Boltzmann equation. For more complete exposure to Boltzmann 
equations see [12] . 

7.3.1 Transport Equation without collisions: 

We have defined particle 4-flow, 



The time component (/i = 0) is the particle density and the space com- 
ponents (/i = 1,2,3) are the particle flow, both measured with respect to 
the observer frame of reference. If dcx^ is a oriented three surface element of 
a plane space-like surface (i.e. a surface whose tangent vector is time) and 
A^a is a small segment situated as x, then we can contract A^^ with da^ to 
obtain the scalar quantity. 




(7.31) 
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ANix) = [ dalN^{x)= [ [d'a,^p'^f{x,p) (7.32) 
In Lorentz frame, d^a^ = {d^x, 0, 0, 0) is purely time like and, 

AN{x) = [ I d^xd^pf{x,p) (7.33) 

Ja'^x J 

AN{x) is just the number of particles in a volume element A'^x In Minkowski 
space, particle is represented by a world line. AA^ thus represent the av. no. 
of world lines crossing the 3-segment A^a. 

ANix,p)= [ [ d'afil-p>^f{x,p) (7.34) 

AN{x,p): the av. no. of particle world lines crossing a segment A'^o" 
with momenta in the range A'^p around p. 

Some time later, same particles will cross a surface element A^a'. Then we 
have the identity, 

/ / d?a,'^p^f{x,p)= j [ d'a,^p'^fix,p) (7.35) 

Ja'^uJa'^P P J A'-^a' J A'-^p P 

Consider the 4-volume A^x enclosed by surface A^a and A^a' and the 
surface of the tube of the world lines. No particle world lines cross the tube 
surface. Thus net flow of particle through the surface A^x of 4-volume A^x 
vanishes. 

/ / d'cr,%'^f{x,p) = 
J A^(T Ja^p P 

Apply Gauss theorem, 

f f d'x^p'^d,f{x,p) = 
Ja^x Ja^p P^ 

Since A^x and A^o" arbitrary 



(7.36) 



(7.37) 



p^d,f{x,p) = (7.38) 
This is the Boltzmann transport equation for a coUisionless system. 
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7.3.2 Transport equation with collisions: 

Number of particles in ranges A^x and A^p changes due to collisions. The 
amount of change can be written as, 

A^x^Cix,p) (7.39) 

where C{x,p) is an invariant function, whose form is to be found. We 
make the following assumptions, 

(a) Only two-particle interactions (dilute system) 

(b) Molecular Chaos (absence of particle correlations). 

(c) f{x,p) vary slowly in space-time. 
Consider a collision of two particles 

According to molecular chaos hypothesis, av. no. of such collisions in 
A^x is proportional to: 

(fpif{x,pi) (fp2 f{x,P2)d!^X 

The proportionality factor, 

W{piP2\p'iP'2) 

is the transition rate. Note we have neglected difference in space-time 
coordinate in f{x,pi) and f{x,p2) (assumption (c) slow variation of f{x,p)). 
The average number of particles lost through collisions is then, 

U^x^ [ ^^^f{x,p,)f{x,p2) X WMp'A) (7.40) 
2 Pi J P2 Pi p'l 

In a similar manner, gain term due to restituting collisions can be calcu- 
lated. 

Boltzmann Equation with collision is 

W = \ j'^^^^^UVlW{p[p',\piP2)-hf2W{pi,P2W^ 
J P2 Pi P 2 

= C[f] (7.41) 
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The transition rate in Eq J7.41l can be related with cross-section of the 
process, 1 + 2 — 2 + 3 and the colhsion term C[f] can be represented in 
several other forms. There is established procedure for solving Boltzmann 
transport equation, e.g. Chapman-Enskog method or Grad's 14-moment 
method. However, in the present course, those methods will not be discussed. 

7.3.3 H-theorem: 

Boltzmann transport equation is manifestly time irreversible. However, the 
microscopic interactions are reversible. How this qualitative change came? 
One of the major assumption in transport equation is the molecular chaos 
hypothesis: colliding particles are un-correlated. This hypothesis makes a 
distinction of past and future. Relativistic transport equation thus describe 
a irreversible process. This property is manifests most clearly in Boltzmann 
H-theorem. In simple terms, H-theorem (or 2nd law of thermodynamics) 
states that entropy production at any space-time point is never negative, 
d^S" > 0. 

H-theorem also defines the equilibrium state of a macroscopic system. In 
the equilibrium state, 

d^S^" = (7.42) 

This condition with transport equation determine the equilibrium distri- 
bution function, it is called Juttner distribution. 

7.3.4 Equilibrium distribution function 

If the macroscopic system is in local thermal equilibrium, at each space- 
time point X, we can specify, in addition to hydrodynamic velocity v{x), a 
temperature T{x) and for each particle species a chemical potential /Xj(x), 
which control the particle density at x. Equilibrium distribution function 
can be obtained from kinetic theory under the condition that df^s^ = 0, i.e. 
it is the distribution which extremises the entropy-four flow. 

Lorentz-covariant equilibrium distribution function can be written as. 



/i,eg('^5 P) 



9i 

g[pM„^(x)-/ii(x)]/r(x) _|_ I 



n+l ~n[p.u(x)- ^i{x)]/T(x} 



(7.43) 
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gi is the degeneracy of the particle, the factor p.u in the exponent is 
the energy of the particle in the local rest frame {p.u p^ = E when 

— )■ (1,0)). The plus and minus sign in the denominator accounts for 
proper quantum statistics of the particle species, (+) for fermions and (— ) 
for bosons. 

7.4 Conservation equations: 

An important property of the collision term in Boltzmann equation is 



with ip{x^p) = a{x) + h^{x)p^. ip{x,p) is generally called summational 
invariant. 

Eq J7.44l can be used to derive conservation equations. For example, con- 
sider the summational invariant il){x,p) = a{x). One obtain. 



Since, we have defined particle 4-current as, A^^ = / ^p'^f{x,p), above 
equation can be written as the macroscopic conservation law of total particle 
number. 



In a system where number of particles of each component is conserved 
separately, one can write. 




(7.44) 








d^N^" = 



(7.45) 



d^Ni: = 0,K = 1,2,. ..N. 
For summational invariant tp{x) = h^{x)p^, one obtain 



(7.46) 
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Energy- momentum tensor is defined as, T'^'^ = J ^p^p'^ f{x) and above 
equation then gives the energy-momentum conservation law, 

d^T^"" = (7.47) 
For a system with a singly conserved charge, the five equations. 



d^N'' = 0, (7.48) 
a^T^^ = (7.49) 

govern the motion of the fiuid. They are called hydrodynamic equations. 
They must be supplemented by the H-theorem or the 2nd law of thermody- 



df,S^ > (7.50) 

Explicit decomposition of energy-momentum tensor and particle four fiow 
are given earlier. For completeness, I repeat them here. 
In Landau-Lifshitz frame, 



= nu^ + ^ (7.51) 
h 

T^"" = eu'^u'' -pA^"" + W (7.52) 

In Eckart's frame. 



iV^ = nu^" (7.53) 

T"" = eu^'u'' -pA^"" + [q^'u^ + q''u''] + W (7.54) 

It is convenient to split the pressure tensor 11'^'" into a traceless part n^'^ 
and the reminder, 

= tt'"^ - nA^'^ (7.55) 

The traceless part tt'^'^ is the shear stress tensor. The reminder 11 is the 
bulk viscous pressure. 
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When the energy-momentum tensor of a closed macroscopic system con- 
tains the dissipative terms, q'^, U, vr^'^, the system is said to be in non- 
equihbrium. The one body distribution function deviates from the equihb- 
rium distribution function given in Eq J7.43[ In the course of time, system 
approaches equihbrium. This approach is dominated by two concepts, ther- 
modynamic forces (the derivatives of macroscopic variables characterising 
the system) and dissipative flows, heat flow {q^), viscous flow U^'^. Flows 
tend to reduce the non-uniformities in the system. Phenomenologically, to a 
good approximation, flows are linearly related to the thermodynamic forces. 
The proportionality constants are called transport coefficients. In ideal or 
inviscid fluid approximation, transport coefficients are exactly zero. Entropy 
is maximised dfj^S^ = 0. Ideal fluid approximations have been widely used 
to model relativistic heavy ion collisions. For ideal fluid, the particle 4-flow 
and energy-momentum tensor can be decomposed as, 

Nf" = nu'' (7.56) 
T^- = (e + p)uV - pg"'' (7.57) 

where, n is the particle density, e is the energy density, p hydrostatic pressure 
and u the four velocity, with the constraint, u"^ = 1.. One notes that the 5 
conservation equations contain 6 unknowns. They are the density n, energy 
density e, pressure p and 3-components of hydrodynamic velocity u { u is 
constrained by definition, = 1). System is closed only with an equation 
of state, p = p{e,n). Equation of state in an important input of hydrody- 
namic calculations. Note that hydrodynamic equations are macroscopic in 
nature. The equations do not depend on the constituents of the system. 
Only the equation of state connects the macroscopic state to the microscopic 
constituents. A hydrodynamic model has the advantage that one can include 
the phenomena of phase transition through the equation of state. 

In general, the fluid is not an ideal. Energy-momentum tensor contains 
the dissipative flows q'^, U and vr'^'^. In kinetic theory, it is possible to relate 
the dissipative flows with gradients of state variables i.e. thermodynamic 
forces. Phenomenologically, if the departure from ideal fluid is small, entropy 
current can be expanded in terms of small deviations, {SN = N'^ — N^, 
ST'^'^ = T'"^ — Tq'^). The subscript denotes the equilibrium value. If the 
expansion contains only the terms flrst order in 6N and 6T'''^ one obtain the 
'flrst order' theory of dissipative hydrodynamics. The procedure yield the 
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'constitutive' relations for heat flow, bulk viscous pressure and shear stress 
tensor, 



n = -CV^n^ 
f = k(V^T- 



T 



TT 



e + p 
27] < Vu" > 



where the angular bracket define a symmetric, traceless tensor. 



< W > 



(7.58) 
(7.59) 
(7.60) 

(7.61) 



In Eqs J7.58l 17.59117.601 ( is the bulk viscosity coefficient, k the heat con- 
ductivity and 7] is the shear viscosity coefficient. The 'constitutive' relations 
can be inserted in the energy-momentum tensor and hydrodynamic equa- 
tions Eqs J7.51|7.52l or Eqs J7.53|7.54l can be solved in principle. However, ffist 
order theory of dissipative hydrodynamics suffers from the problem of causal- 
ity. They are acausal. Also, there may be the problem of instability. The 
causality problem is removed if the expansion of entropy current contains 
terms in second order in deviations 6N and ST'^'^. The procedure gives rise 
to the relaxation equations for the dissipative flows. In addition to the con- 
servation equations, relaxation equations for the dissipative flows need to be 
solved. 



n = -CVu^-mDU (7.62) 

gA* = k(V^T —V^p)-r„Dq'' (7.63) 

e + p 

n"" = 27]< W^'u" > -T^Dn^"" (7.64) 

D = u^d^ is the convective time derivative, rn, Tq and t^^ are relax- 
ation time for the bulk viscous pressure, conductivity and shear stress tensor. 
Transport coefficients and as well as the relaxation times, in principle, can 
be calculated in kinetic theory. However, for strongly interacting system like 
QGP, calculations are very complex. I may note here that I have written 
the simplest possible relaxation equation. Relaxation equations can contain 
additional terms. 
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There are 10 independent dissipative flows. The shear stress tensor tt^" is 
symmetric, traceless and transverse to hydrodynamic velocity (n^jTr^*^ = 0). 
TT^'^ has 6 independent components. Heat flow is also transverse to hydro- 
dynamic velocity {u^w^ = 0) and has 3 independent components. And there 
is the bulk viscous pressure. In second-order hydrodynamics, in addition to 
the 5 conservation equations, 10 relaxation equations for dissipative flows 
need to be solved. 

Above discussions about dissipative hydrodynamics are rather sketchy. 
For more information one can see [17] [l8] [l9] [50] [51] [52] [53] . 

8 Hydrodynamic model for heavy ion colli- 
sions 

8.1 Different stages of HI collisions: 

We want to study properties of QGP. QGP existed in early universe. It is 
also possible that it exists at the core of neutron stars. However, QGP at 
the early universe or at the core of neutron stars are not accessible for study. 
Theoretical considerations led us to believe that QGP can be produced in 
laboratory by colliding heavy nuclei at very high energy. A nucleus-nucleus 
collision is a well established tool to study properties of nuclear matter. 
Recent experiments at Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory and Large Hadron Collider (LHC) at CERN, strongly 
suggests that QGP is formed in high energy nuclear collisions. 

Let us qualitatively discuss the collision process with increasing collision 
energy. In very low energy collisions, nucleus as a whole interacts. Indication 
are obtained from, say giant dipole excitations, where the compound nucleus 
undergoes dipole oscillation (proton and neutron fluid oscillate out of phase) . 
One can also excite the nucleus and populate various excited states. As the 
energy is increased, nucleons in the nucleus start to interact, one can see 
production of new particles, e.g. vr. At still higher energy, quarks inside the 
nucleons will interact. Here also production of different particle species will 
be observed. However, in contrast to low/medium energy nuclear reactions, 
where one can describe pA/AA collisions entirely in terms of NN collisions, 
in relativistic energy, such a description will fail. 

A nucleus-nucleus collision at relativistic energy passes through different 
stages. Initial collisions are expected to be in the partonic level. Schematic 
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Figure 26: different stages of nuclear collisions. 

picture of different staged of the collisions are shown in Figl26l One can 
broadly classify different stages; 

(i) Pre-equilibrium stage: initial partonic collisions produce a fireball in 
a highly excited state. In all possibility, the fireball is not in equilibrium. 
Constituents of the system collide frequently to establish a 'local' equilibrium' 
state. The time takes to establish local equilibrium is called thermalisation 
time. It is an important parameter. Explicit comparison of hydrodynamic 
models with experiments indicate that in A/ijy^=200 GeV Au+Au collisions, 
local equilibrium is established sufficiently fast, in the time scale 0.5-1.0 fm. 
However, it is not understood how such fast equilibration can be achieved. 
In classical electromagnetic plasma, self excited transverse modes exist in 
plasmas with momentum anisotropy [5l]. They are called Weibel instability. 
The instability can grow very fast. The instability reduces the momentum 
anisotropy. It is believed that in QGP, non-abelian version of the Weibel 
instability establishes rapid thermalisation [55] . 

(ii) expansion stage: In the equilibrium or the thermalised state, the 
system has thermal pressure which acts against the surrounding vacuum. The 
system then undergoes collective (hydrodynamic) expansion. As the system 
expands, its density (energy density) decreases and the system cool. The 
expansion and cooling is governed by the energy-momentum conservation 
equations, which involve an equation of state p = p{e,nB)- Now if there is 
a QCD phase transition in the model, then below the critical energy density 
ecr ~ IGeV/ fm^, or critical temperature Tcr ~ 200MeV, the partons (quarks 
and gluons) will convert to hadrons. In the hadronisation stage, over a small 
temperature interval, entropy density will decrease very fast. Since total 
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entropy can not decrease, it implies that the fire ball will expand rapidly, 
while temperature remains approximately constant. If the transition is 1st 
order, there will be mixed phase, where QGP and hadronic resonance gas can 
co-exist. In the mixed phase, speed of sound = ^ ~ 0. In 2nd or cross-over 
transition, there will not be any mixed phase, even then, near the transition, 
equation of state is soft and speed of sound is small (see Fig JT7|) . Then 
collective flow will not grow much during the transition period. Ultimately 
all the partonic matter will be converted into hadronic matter. 

(iii) Freeze-out: Hadronic matter will also be in thermal equilibrium. 
Constituent hadrons will collide to maintain local equilibrium. The system 
will expand and cool. A stage will come when inelastic collisions, in which 
hadron changes identity, become too small to keep up with expansion. The 
stage is called chemical freeze-out. Hadron abundances will remain fixed after 
the chemical freeze-out. However, due to elastic collisions, local equilibrium 
can still be maintained and system will cool and expands with fixed hadron 
abundances. Eventually a stage will comes when average distance between 
the constituents will be larger than strong interaction range. Collisions be- 
tween the constituents will be so infrequent that 'local' thermal equilibrium 
can not be maintained. The hydrodynamic description will break down. The 
hadrons decouple or freeze-out. It is called kinetic freeze-out. Hadrons from 
the freeze-out surface will be detected in the detector. 

Hydrodynamics provides a simple, intuitive description of relativistic 
heavy-ion collisions. Hydrodynamic models requires the assumption of lo- 
cal (thermal) equilibrium, i.e. one assume that at each space time point x 
of the fluid, a small region can be considered where equilibrium is achieved, 
a temperature T(x) can be defined. As discussed above, from the equilib- 
rium stage to the kinetic freeze-out, relativistic heavy ion collisions can be 
modeled by hydrodynamics. 

A variety of experimental data from Relativistic Heavy Ion Collider (RHIC) 
experiments, are successfully explained in ideal hydrodynamical models [56] [57] 
Some problems however remained. For example, ideal hydrodynamics de- 
scription of experimental data becomes poorer as the collisions become more 
and more peripheral. Also, contrary to the experimental data, in ideal hy- 
drodynamics, elliptic flow (it is an important observable in high energy heavy 
ion collisions and will be discussed later) continue to increase with transverse 
momentum. Deficiencies of ideal hydrodynamic models are to some extent 
corrected in dissipative hydrodynamics [58] [59] . 
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8.2 Hydrodynamic equations 

In this section, I describe the general procedure followed in hydrodynamic 
modeling of heavy ion collisions. For simplicity, we assume an ideal fluid 
with a single conserved charge (e.g. baryon density). The five conservation 
equations. 



9^iV^ = 0, (8.1) 
= (8.2) 

govern the motion of the fluid. Given an initial configuration of the fluid and 
an equation of state p = p{e, ub), the equations can be solved numerically to 
obtain the space-time evolution of the fluid. 

In heavy ion collisions, appropriate coordinates for solving hydrodynamic 
equations are (r, x, y, rj) rather than (t, x, y, z). 



x" = {t,x,y,z) x"" = (r, x, y, r/) (8.3) 
t = rcosh?7; r = \Jt'^ — z'^ (8-4) 



z = rsinh?7 : 




(8.5) 



In (r, X, y, rf) coordinate system the metric is. 



ds^ = g^vdx^dx" = dt^ — dx^ — dy"^ — dz'^ (8.6) 

= dr"^ - dx^ - dy^ - T^dri'^ (8.7) 

and , 

g^'' = dzag{l,-l,-l,-l/T^). (8.8) 

One note that the space-time is not flat anymore, it is curved. Accord- 
ingly, one need 'affine connections' or the Christoffel symbols, 

2^ \dx'' ^ dx' dx^J ^ ^ 
In {T,x,j,'r]) coordinate only non-zero Christoffel symbols are. 
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r;, = r,n^ = l/r (8.10) 
Covariant (semicolon) derivative of a contravariant tensor is given by, 



Ai — i_ r« Am 

f)Aik 

Aik I pi Amk I -pk Ami 

Five conservation Equations in {T,x,j,'r]) coordinate system can be easily 
derived. For 'ideal fluid', they are. 



= {drN^ + r;,„Ar™) + 9,iv" + dyNy + [dr^N'^ + r;|„Ar™ 

= drN^ + +d^N'' + dyNy + d^N'' + -N^ (8.11) 
(n)T^- = = + + + 

= (9,T"" + 2r;^T'"") + a^.T™ + dyry 

= 9,T"" + 9^.T™ + (9yT"^ + S^T"" + rT"^ + ^T"" (8.12) 
(m)T.^- = = + T.^" + T.7 + T.:j" 

= d^T^y + d^T''y + dyTyy + dr,Ty'^ + -T^'' (8.13) 
(it;)T.;i?' = = T^y + T^y + T^y + T;^y 

= drT^y + d^T'^y + dyjyy + dr^T^^ + -T^y (8.14) 
(t;)T.7 = o = T^;^ + T^j' + Ty'^ + t;i;' 

= drT^"^ + d^T'''^ + dyTy + d^T'T^ + ^T'^^ (8.15) 

In ideal hydrodynamics, or in 1st order hydrodynamics with dissipation, 
the five partial equations noted above are to be solved simultaneously to ob- 
tain space-time evolution of the fluid. Flux corrected SHASTA algorithm [60] 
can be used to solve those equations. In 2nd order dissipative hydrodynamics 
however, 10 additional relaxation equations have to be solved simultaneously. 
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For illustrative purpose, let us specialize in one-dimensional, Bjorken scal- 
ing flow [61]. In one dimension expansion, hydrodynamic four velocity is 
= (1,0,0,0) and relevant energy-momentum components are. 

Inserting the values in Eq J8.12l we obtain, 

= drT^^ + d^T^"" + dyT^y + drfT^'" + tT"" + ^T^^ 

= drT^^ + rT'"' + -T"" = ^ + (8.17) 

T OT T 

Using entropy density s = the Eq J8.17l can be recast in to 

(Its 



0, (8.1^ 



dr 

which can be solved as, st = constant. One dimensional flow is isentropic. 
In ideal gas, s oc and we get the well known law for Bjorken scaling 
expansion, 

TS = TfTf (8.19) 

Given the fluid temperature at initial time r^, its value at a subsequent 
time is easily obtained. Similarly, one can also solve for the number conser- 
vation equation. For fluid velocity u'^ = (1,0,0,0), A^'^ = (n, 0,0,0). From 
Eq J8.11l we obtain, 

= drN^ + +d^N'' + dyNy + dr^N"^ + ^N^ 

= drn+-n (8.20) 
r 

The Eq. can be solved to give, 

nf = n— (8.21) 

Noting that density is inversely proportional to volume, one find the scal- 
ing law for volume expansion, 

Vf = V,^ (8.22) 
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8.3 Cooper-Frye prescription for particle distribution 

Hydrodynamic equations give the space-time evolution of the fluid till a 
freeze-out condition. One obtain information about the fluid energy den- 
sity (or temperature) and velocity at the freeze-out. The information needs 
to be converted into particle distribution. This is done through Cooper-Frye 
prescription [62] . 

Consider a three dimensional hypersurface S(x) in 4-dimensional space- 
time and count the number of particles crossing the hypersurface. Let da^ is 
an infinitesimal element perpendicular to S(x) and directing outward. If 
is the current of particles, then the scalar product da^j^, gives the number 
of particles crossing the infinitesimal surface da. The total number crossing 
the hyper surface S is, 



N 



da,f = ^ da, I ^p'^f{x,p)^ (8.23) 



with f{x,p) the one body distribution function. 
In differential form. 

This is the Cooper-Frye prescription [62] for the invariant distribution of 
particles. Let us evaluate the term da^p^. 

In (r, X, y, rj) coordinates, the freeze-out surface can be parameterised as, 

cr'^ = (r/(x, y) cosh?7, x, y, Tf(x, y) sinhr/), (8.25) 
and the normal vector on the hyper surface is, 

OTf dTf 

da'^ = (coshr?, — -^r^, — smh.r])Tfdxdydr], (8.26) 
ox ay 

4-momentum p^ can be parameterised as, 

p^ = {rriT cosh y, Px,Py, triT sinh y) (8.27) 
The volume element 

dTf dTf 

p^da^, = {rriT cosh(?/ -7])- Px-^^ - Py-^)Tfdxdydri (8.28) 
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/ in 18.241 is the distribution function. In ideal hydrodynamics, the / 
corresponds to Lorentz-covariant equihbrium distribution function, 



9i 



^[pi^u^{x)-iii{x)]/T{x) I 
oo 

= ^. ^(^ip-)"+lg-"[p"«M(^)-M{a:)]/T{a;) (8.29) 
i=l 

Qi is the degeneracy of the particle, the factor p.u in the exponent is 
the energy of the particle in the local rest frame {p^u^ = E when 

— )■ (1,0)). The plus and minus sign in the denominator accounts for 
proper quantum statistics of the particle species, (+) for fermions and (— ) 
for bosons. 

Let us evaluate the term p^u^. The fluid 4-velocity can be parameterised 

as, 



u 



= 7r(cosh?7, Wj;, Wj^, sinh?7), (8.30) 



VI - 4 V 1 - 4 - 

The scalar product p.u in the equilibrium distribution function is then 
easily calculated as, 

P^u^, = 'jrimT cosh(?/ - rj) - v^p^ - VyPy) (8.32) 

Eqs J8.28|8.32l completely specifies the Cooper- Frye invariant distribution 
EqJOl 

In dissipative hydrodynamics, distribution function / deviates from the 
equilibrium value, 

/(x,p) = /ejl + 0(x,p)] (8.33) 

Since distribution function is a scalar, the deviation function must be 
written as sum of scalars constructed from p^, and dissipative flows. It 
can be parameterised as, 

0(a;, p) = A{x,p)Ii + B{x,p)p^,q^ + C{x, p)p^p^n^"' (8.34) 
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The coefficients A, B and C can be determined from the condition of 
the fit, that the number density and energy density is determined by the 
equihbrium distribution functions, 



^p%/e,0 = (8.35) 



/ ^(p"«.)Ve,0 = (8.36) 

I will not discuss in detail, but the deviation function for shear viscosity 
can be written as, 



(j^shear = C {x , p)p^p^7rf"' = ^ ^ ^ ^^^^ P^P^.TT^" (8.37) 

In [63], deviation function for bulk viscosity is obtained from Grads 14 
moment method. 



= A{x, p)U = Dopf^u^ + Bop'^p'^A^, + Bop'^p'^u^u, (8.38) 
expressions for Dq, Bq and Bq can be found in [63] . 



8.4 Initial conditions for hydrodynamic analysis 

One understands that hydrodynamics is an initial value problem. For exam- 
ple, in Bjorken one dimensional hydrodynamics, given the density/temperature 
at some initial time Tj (the time beyond which hydrodynamics is applicable, 
itself a parameter of the model), density/temperature evolution of the fluid 
can be obtained. A kinetic freeze-out condition also required to define the 
freeze-out surface such that Cooper-Frye prescription can yield particle's in- 
variant distribution. One simple procedure to implement kinetic freeze-out 
is to assume a fixed freeze-out temperature Tp. Depending on the model, 
experimental data are fitted with Tp=100-140 MeV. 

In a more general system, one has to initialise the (baryon) number den- 
sity n{x,y,ri), energy density e{x,y,r]) and velocity u{x,y,r]) = -y\{x,y,ri) 
distributions at the initial time r^. Indeed, one of the aims of hydrodynamic 
analysis of heavy ion collisions in ultra-relativistic collisions is to obtain the 
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initial conditions of the produced fluid, by comparing hydrodynamic simu- 
lations with experimental data. As discussed earlier, experimental results 
are given in terms of collision centrality. One then tries to parameterise the 
initial condition in terms of impact parameter, such that once the parame- 
ters are flxed at some particular collision centrality, it can predict for other 
collision centralities. 

In (r,x,y,?7) coordinate, for the initial energy density, a common practice 
is to assume a factorised form. 



B{x,y) being the initial energy density in the transverse plane and -^(77) in 
the direction of (spatial) rapidity. One can use a Gaussian distribution for 
H[r)). Transverse energy distribution e{x,y) can be conveniently parame- 
terised in a Glauber model or in color glass condensate (CGC) model. The 
number density distribution or the velocity distribution at the initial time can 
be similarly parameterised. In general, one assume zero initial fluid velocity 
at the initial time, though it is possible that fluid have non-zero velocity, 
especially near the surface of the fluid. The reasoning is simple. Fluid con- 
stituents can have random velocity. In the interior of the fluid, the random 
velocities will balance to produce net zero velocity. However, near the surface 
random velocities will not be balanced. 

8.4.1 Glauber model initial condition 

In sectionjSl I have discussed the Glauber model. Expressions for the number 
of participant nucleons and number of binary collisions, in impact parameter 
b collisions were obtained. 



From the above equations, , transverse proflle of binary collision number 
and participant numbers, in impact parameter b collisions can be easily 
obtained as. 






(8.40) 



(8.41) 
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N,ou{x,y) = ABaNNTA{x + h/2,y)TB{x-h/2,y) (8.42) 
Npart{x,y) = ATA{x + h/2,y){l-[l-aNNTB{x-h/2,y)f) 

+ BTsix - b/2, y){l - [1 - a^NTAix + b/2, y)f) (8.43) 

Comparison of hydrodynamic simulations with experimental data indi- 
cate that a combined profile, 



e{x,y) oc [il-f)Npart{x,y)+fN^ou{x,y)] = eo[{l-f)Npart{x,y)+fNcoii{x,y)] 

(8.44) 

with / ^ 0.1 — 0.2 best explains the data. Once the proportionality factor 
Co is fixed in a given collision centrality, the impact parameter dependence 
of the model allow one to predict for the energy density distribution at other 
collision centralities. 



8.4.2 CGC initial condition 

In sectionjni I have discussed Color Glass Condensate (CGC). CGC is a 
quantum mechanical state of matter at high energy. The earliest state of 
matter produced in high energy nucleus-nucleus collisions, may not be much 
different from this quantum mechanical state. At very early time, CGC 
evolves in to a distribution of gluons. Later these gluons thermalise and form 
QGP. CGC models have been used extensively to model the transverse energy 
distribution of the initial QGP fluid in hydrodynamic models |M] [53] [HS] In 
the following, I briefly describe the procedure to obtain the initial condition 
in high energy nuclear collisions in the KLN (Kharzeev-Levin-Nardi) [71] 
approach to CGC. 

The number of gluons produced in the factorisation formula is given 

by, 



dNg Att^N^ f (PpT 



(fkTOisikT) 



(PrrdY N"^ - I J Pt 

X (t)A{xi,PT,rT)(t)B{x2,{VT-kTf,rT) (8.45) 

where px and Y are the transverse momentum and rapidity of the produced 
gluon. Xi2 = ^e^^ is the momentum fraction of colliding gluons ladders at 
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cm. energy y/s. as{kT) is the strong coupling constant at the momentum 
scale kx- The unintegrated gluon distribution functions (pA in Eq J8.45l are 
related to the gluon density in a nucleus at the transverse position rr(= x, y), 

xG{x,Q^)= dPkTd'^rT(l)A{x,kl-rT) (8.46) 

In principle, unintegrated gluon distribution function should be a so- 
lution of non-linear quantum evolution equation e.g. JIMWLK equation 
[67] [68] [69] [To] . In Kharzeev-Levin-Nardi (KLN) approach [71] approach 
(which captures the essential features of the gluon saturation), the unin- 
tegrated gluon distribution functions are taken as [7T|[72] . 

where Qs is saturation momentum at the given momentum fraction x and at 
the transverse position tt- 

In the KLN approach, the saturation scale in AB collision is parameterised 

as mm. 




The form Qs{x) ~ x with A ~ 0.2 — 0.3 is motivated from DIS experi- 
ments. Np}^^ in the above equation is the transverse density of participant 
nucleons, which can be calculated in a Glauber model (e.g. see Eq J8.43p . 



N^artirr) = ATa{x + b/2, y)(l - [1 - ar^NTsix - b/2, y)] 



(8.49) 



In the CGC model, the transverse energy density should follow Eq J8.45[ 
However, Eq J8.45l is valid in the time scale ~ when the medium may 
not be in thermal equilibrium. One assumes that the medium undergoes one 
dimensional Bjorken (longitudinal, isentropic) expansion during the period 
Ts to Tj. The density at the time Tj, when hydrodynamic become applicable 
is easily obtained as n(rj) = ^n{Ts). The transverse energy density profile 
at the initial time r,- is then, 



Ax,y,h) = eo 



dxdydY 



n4/3 



^.50) 
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Figure 27: Schematic picture of a non-zero impact parameter collision, 
with Co a normalising factor, which is to be fixed from experimental data. 
8.5 Collective flow 

In relativistic heavy ion collisions, one of the important observables is the 
azimuthal distribution of produced particles. In Figj271 geometry of a col- 
lision at non-zero impact parameter collision is shown. The overlap region 
of the two nuclei is the participant region, where most of the collisions oc- 
cur. The target and projectile remnants on the periphery acts as spectator. 
It is obvious from Figj27t that in non-zero impact parameter collision, the 
participant or the reaction zone in coordinate space do not posses azimuthal 
symmetry. Multiple collisions among the constituent particles translate this 
spatial anisotropy into momentum anisotropy of the produced particles. The 
observed momentum anisotropy is called collective flow and has a natural 
explanation in a hydrodynamic model [57]. In the following we briefly dis- 
cuss collective flow phenomena. More detailed expositions can be found in 
[73] [71] [75]. 

Momentum anisotropy is best studied by decomposing the invariant dis- 
tribution in a Fourier series. For example, the momentum integrated invari- 
ant distribution of a particle can be expanded as, 

^ = ^[1 + 2 Vi;„cos[n(0-?A)],n = 1,2,3... (8.51) 

n 

(j) is the azimuthal angle of the detected particle and i/j is the plane of the 
symmetry of initial collision zone. For smooth initial matter distribution, 
plane of symmetry of the collision zone coincides with the reaction plane 
"^RP (the plane containing the impact parameter and the beam axis). The 
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Figure 28: temporal evolution spatial eccentricity (ex) and momentum 
anisotropy (e^) with time [57]. 

sine terms are not present in the expansion due to symmetry with respect to 
the reaction plane. 

Flow coefficients Vn are easily obtained, 



/ d(l)^cos{n(j) — nip) 

7W 



{cos{n(i) - ni))) = '^^ ^ ^ ,n = 1,2,3... (8.52) 



vi is called (integrated) directed flow, V2 (integrated) elliptic flow, 
(integrated) triangular flow, (integrated) hexadecapole flow etc. 

Similar to Eq J8.51l one can Fourier expand the invariant distribution 

dP'N 1 (PN 

^ W = 2^?^|l + 2 ^ - n = 1, 2, 3... (8.53) 

and differential flow coefficients are obtained as, 



f dd) — f^,,j cosind) — nip) 

J ^'i^ pj,dprpd(f)dy 



VniPT, y) = {cos{n(f) - mjj)) = P^'iPTdHy n = 1, 2, 3... 



(8.54) 

Second flow coefficient has been studied extensively in RHIC and LHC 
energy collisions. Finite-non-zero value of V2 is thought to be direct signature 
of production of thermalised medium. 
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Elliptic flow in heavy ion collisions is best understood in a hydro dynamic 
model [57j. Elliptic flow measure the momentum anisotropy. In non-zero 
impact parameter collisions, the reaction zone is spatially asymmetric. Spa- 
tial asymmetry of the initial reaction zone can be quantifled in terms of 
eccentricity, defined as, 

where (..) indicate energy /entropy density weighted averaging. In non-zero 
impact parameter collision, initial eccentricity is non-zero, positive. If a 
thermalised medium is produced in the reaction zone, due to thermodynamic 
pressure, the medium will expand against the outside vacuum. One can 
immediately see that pressure gradient will be more along the minor axis 
than along the major axis. Due to differential pressure gradient, as the 
system evolves with time, eccentricity will reduce. Momentum distribution 
of particles are isotropic initially. If momentum anisotropy is measured as, 

e,= {fflr-r! (8.66) 

initially Sp will be zero. However, as the fluid evolves, rescattering of parti- 
cles will introduce asymmetry and Sp will grow. It is expected to saturate 
beyond certain time, when reaction zone attains azimuthal symmetry. In 
that sense, elliptic flow is self quenching phenomena, driving force of the 
flow (the reaction zone asymmetry) continuously reduces as the flow grow. 
In Figj28| ideal hydrodynamic model simulations for temporal evolution of 
spatial eccentricity and momentum anisotropy are shown. They follow our 
expectations. 

The second harmonic coefficient or the elliptic flow {V2) has been studied 
extensively in A/i^jy=200 GeV Au+Au collisions at RHIC [SIE]. Recently, 
ALICE collaboration measured elliptic flow in ^/s^^=2.7Q TeV Pb-|-Pb col- 
lisions at LHC [76l[77]. Large elliptic flow has provided compelling evidence 
that at RHIC and LHC, nearly perfect fluid is produced. Deviation from the 
ideal fluid behavior is controlled by shear viscosity to entropy ratio {rj/s). 
Effect of shear viscosity is to dampen the flow coefficients. Elliptic flow 
has sensitive dependence on rj/s. In smooth hydrodynamics, sensitivity of 
elliptic flow has been utilised to obtain phenomenological estimates oi rj/s 
I65l[78l[79l|80l[8ll|82l[83l|8l]. It appears that QGP viscosity over entropy 
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Figure 29: (left panel) Charged particles pseudo rapidity density in 
1/5^^=130 GeV Au+Au collisions, in a central 0-6% and peripheral 35-45% 
collision centralities are compared with hydrodynamic simulation with two 
initial condition A and B (for details of the initial conditions see [86]). (mid- 
dle panel) transverse momentum spectra of identified particles in -y/s^^=130 
GeV Au+Au collisions in different centrality ranges of collisions (see |57] for 
details) and (right panel) differential elliptic flow of identified particles is in 
A/ijy^=200 GeV Au+Au collisions (see [87] for details). 

ratio is close to 77/5 ~ l/Air. 

8.6 Some results of Hydrodynamic simulation of heavy 
ion collisions 

Various authors have simulated heavy ion collisions at relativistic energy. In 
the following I will show some representative results. In the left panel of 
Figl29l PHOBOS measurements [85] of charged particles rapidity density in 
^/sJ^Jy=200 GeV Au+Au collisions, in two different collision centralities are 
shown. Experimental data are nicely reproduced in hydrodynamic simula- 
tions. For details see [86] . 

In the middle panel of Fig|29| fits obtained to the identified particles spec- 
tra in STAR and PHENIX experiments in a hydrodynamic model simulation 
are shown. Note the quality of fit. Data in the collision centrality 0-5% to 
60-90% are well explained. Finally, in the right panel of Figj29| elliptic flow 
of identified particles, as a function of pr are shown. One observe that exper- 
imental flow of identified particles are mass dependent, more flow for lighter 
particles than for heavier particles. It is called mass splitting of flow. At 
high pt however, effect of mass splitting is reduced. In Figj29l the solid lines 
are hydrodynamic model simulations. Experimental mass splitting of flow 
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Table 6: Central energy density [si) and temperature (Tj) at the initial time 
Ti=O.Q fm/c, for different values of viscosity to entropy ratio (rj/s). The 
bracketed values are estimated central energy density and temperature in 
^/s^Jy=200 GeV Au+Au collisions [79]. Initial time of the simulations is 
Ti=0.6 fm. 



r]/s 





0.08 


0.12 


0.16 




89.2 ±5.0 


78.0 ±4.0 


70.5 ±3.5 


61.7±3.0 




(35.5 ±5.0) 


(29.1 ±3.6) 


(25.6 ±4.0) 


(20.8 ±2.7) 




486 ±6 


475 ±5 


462 ±6 


447 ±5 


(MeV) 


(377 ± 14) 


(359 ± 12) 


(348 ± 14) 


(331 ± 11) 



is correctly reproduced in hydrodynamic simulations. The agreement with 
experiment is also good. The simulations results shown in FigJ^Hlare for ideal 
fluid only. As mentioned earlier, there are several simulations with viscous 
fluid. I will not show the results. I just mention that compared to ideal fluid, 
a viscous fluid will require less initial energy density or temperature. This 
is because, entropy is generated during viscous evolution. As an example, 
in tablejni I have noted the central energy density and temperature of the 
fluid obtained from fits to experimental data in A/i^^=200 GeV Au±Au 
and y/s^j^=2.7Q TeV Pb±Pb collisions. Note that viscous fluid require less 
energy density or temperature. See [79], [81] for details of the simulations. 

8.7 Event by event hydrodynamics 

Very large number of particles in the final state, in RHIC and LHC energy 
collisions, enables experimentalist to analyse the experimental data event- 
by-event. Importance of event-by-event analysis is best explained by the 
classic example by A. D. Jackson. If in a rainy day, you hold out a sheet 
of paper outside the window and forget about it for a long time, you find 
it uniformly soaked. You conclude that the spatial distribution of rain was 
uniform. However, if you continue to look into the paper, you find that spatial 
distribution is not uniform in a short time scale. Similarly, high statistics 
data in a single event may be very different from the data averaged over many 
events and reveal interesting physics. More details about event-by-event can 
be found in [88] . 

One of the aims of hydrodynamic analysis is to find the initial conditions 
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Figure 30: Schematic picture of Monte Carlo Glauber simulation for partici- 
pating nucleons in the transverse plane. The overlap region of participating 
nucleons are tilted with respect to the reaction plane. 

of the matter produced in heavy ion collisions. Now the initial conditions can 
fluctuate, event by event. In Figj30l a schematic picture for Glauber Monte- 
Carlo simulation for the participating nucleons in the transverse plane is 
shown. The overlap region of participating nucleons is tilted with respect 
to the reaction plane. It is an important realisation that azimuth of the 
particles should be measured with respect to the participating plane angle 
rather than the reaction plane angle. 

Now, the positions of the participating nucleons will fluctuate event-be- 
event, so does the participating plane angle. The participating nucleons 
which determine the symmetry plane (ippp), will fluctuate around the re- 
action plane [89]. As a result odd harmonics, which were exactly zero for 
smoothed initial distribution, can be developed. Recently, ALICE collabora- 
tion has observed odd harmonic flows in A/i^^=2.76 TeV Pb-|-Pb collisions 
[77] . In most central collisions, the elliptic flow (^2) and triangular flow (v-^) 
are of similar magnitude. In peripheral collisions however, elliptic flow domi- 
nates. Triangular flow is being investigated in event-by-event hydrodynamics 
[90], [91], [92], [93], [91], [95], [96], [97]. It has been conjectured that third har- 
monic W3, which is response of the initial triangularity of the medium, is 
responsible for the observed structures in two particle correlation in Au+Au 
collisions. The ridge structure in pp collisions also has a natural explanation 
if odd harmonic flow develops. It is also expected that triangular flow will be 
more sensitive to dissipative effects and can constrain the phenomenological 
estimate of shear viscosity to entropy ratio. However, detailed simulations in- 
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dicate that in event-by-event hydrodynamics, sensitivity of V2 and ^3 toward 
viscosity is reduced due to fluctuations [96].[97]. 

9 Signals of Quark-Gluon-Plasma 

We expect that QGP is produced in uhra-relativistic heavy ion colhsions. 
However, how do we know it is produced indeed, or, which observable will 
signal formation of QGP. Signals of QGP are very important subject in the 
study of QGP. Unlike in other phase transitions, in conflnement-deconflnement 
phase transition (or cross over), the constituents of the high temperature 
phase (QGP) do not exist freely. They are conflned within the hadrons. The 
problem is closely related to quark conflnement; quarks are unobservable. 
QGP, even if produced in a collision, is a transient state, it expands, cools, 
hadronises, cools further till interactions between the hadrons become too 
weak to continue the evolution. Any information about the QGP phase, if 
produced in high energy nuclear collisions, has to be obtained from the ob- 
served hadrons only. Hadronisation is a non-perturbative process. Till now, 
it is not properly understood. Whether or not the hadronisation process 
erases any memory of the constituent quarks is uncertain. If the hadroni- 
sation process erases the memory, from the observed hadrons one can not 
comment on the initial QGP phase. Present search for QGP at RHIC/LHC 
is on the premise that the hadronisation process does not erase the memory 
and from the observed hadrons, one can comment on the possible existence 
of QGP. In the following, I briefly discuss some important QGP signals. 

Table 7: J/'0 mass, radii and formation time from solution of Schrodinger 
Equation. 





J/ij 




Xc 


M(GeV) 


3.07 


3.698 


3.5 


R(rm) 


0.453 


0.875 


0.696 


7-F(fm) 


0.89 


1.5 


2.0 



9.1 J/i/j suppression 

Charmoniums (bottoniums) are bound state of cc (bb) quarks. Charm and 
bottom quarks are heavy {rricharm ~ 1.15-1.35 GeV, rrihottom ~ 4.0-4.0 GeV) 
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Figure 31: quark potential as a function of temperature. 



and non-relativistic Schrodinger equation can be solved to obtain the bound 
state properties, with inter-quark potential (Cornell potential) 



V{r) 



a 



ar 



eff 



(9.1) 



where r is the inter quark separation, a ~ 0.192 GeV^ is the string constant 
and a^ff ~0.471. Quarkonium mass, radius and formation time can be ob- 
tained from solution of non-relativistic Schrodinger equation and are given in 
the tableJTl Now at high temperature, interaction potential will be screened, 



V{r,T) 



a 



/i(T) 



(1 _ e-^(^)'-) 



r 



(9.2) 



/i(T) is the inverse of the screening radius (Debye radius) and is called the 
screening mass. For /i — )■ 0, Eq J9.1l is recovered. For /i 7^ 0, the screened 
potential satisfies. 



Limr^Q[rV{r,T)] ~ —a 
the 1/r behavior in the short distance limit. For large r. 



Lirrir 



-In 



a 



-V{r,T) 



(9.3) 



(9.4) 



indicating that the range of the binding force decreases exponentially with 
screening mass. In Figj3T| inter-quark potential as a function of temperature 
is shown. The screening mass /i is as increasing function of temperature T. 



From lattice simulations, ^ A^. 

J-c J-c 



In 1986, Matsui and Satz [HH] suggested that if QGP is formed in nuclear 
collisions, J/?/', the bound state of cc will be suppressed, w.r.t. pp collisions. 
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Figure 32: Nuclear modification factor Raa for J/tp in A/i^^=200 GeV 
Au+Au collisions. 

The idea is simple. In presence of QGP, J/ip production will be inhibited 
due to screening of potential. A cc pair, which could transform into a, J/ip is 
now unable to do so. Over the years, several groups have measured the J/ip 
yield in heavy ion collisions (for a review of the data prior to RHIC energy 
collisions, and the interpretations see Refs. [991 llOOj ). In brief, experimental 
data do show suppression. However, suppression is observed in pA collisions 
also, where, one does not expect QGP formation. It is understood that in an 
inelastic collision with nucleons, J/ip's can be dissociated and lead to sup- 
pression. Suppression in pA collisions is termed cold nuclear matter (CNM) 
effect. It is important to disentangle CNM effect from the experimental data 
to obtain the suppression due to deconfinement. 

PHENIX collaboration has made systematic measurements of J/ip pro- 
duction in nuclear collisions. They have measured J/ijj yield in p+p collisions 
at RHIC to obtain the basic 'reference' invariant yield |10H 1102"] . Measure- 
ments of J/ip production in d+Au collisions |102j give reference for cold 
nuclear matter effects. Cold and hot nuclear matter effects are studied in 
Au+Au and Cu+Cu collisions, where yields are measured as a function of 
collision centrality |103[ 1104] . In Figj32| PHENIX measurements for nuclear 
modification factor. 



98 




Figure 33: STAR measurements for enhancement factor in y/sj^j^=200 GeV 
Cu+Cu and Au+Au collisions, for various multistrange particles. 



I ^AA^J/i(,X 

in -y/ijy^=200 GeV Au+Au collisions are shown. Data shows suppression. 
Only CNM effect can not explain the data. If J/ip^s are suppressed in QGP, 
data are explained. At RHIC energy, it has been argued that rather than 
suppression, charmonium's will be enhanced |105[ 1106] . Due to large initial 
energy, large number of cc pairs will be produced in initial hard scatterings. 
Recombination of cc can occur enhancing the charmonium production. Ap- 
parently, PHENIX data on J ftp production in Au+Au are not consistent 
with models which predict J/ip enhancement. 



9.2 Strangeness enhancement 

For long, strangeness enhancement is considered cLS cL SI gnature of QGP for- 
mation |1U7] . In QGP, strangeness will be more abundant that non-strange 
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quarks. For example, in equilibrating plasma, strange quark density (with 
degeneracy 2{spin) x 3{color), can be obtained as (see Eq J4.14p . 



= = 3 X 2 y ^e-v^+"^^/^ = ^^K,{mjT) (9.6) 

density of non-strange anti-quarks q {q stands for u or d, in the limit of 
small mass, is 

^_ = 3 X 2 1 -0-^e-\^\/^e-^^/^ = e'^^/^ (9.7) 

where /Xg = is the quark chemical potential. The ratio of strange quarks 
to non-strange quarks, 

^ = im'A,me«/- (9.8) 



3Tm 



2 



riq 2\T J \T 

For /ifi > 0, the ratio is greater than one. Strange quarks will produce in 
more abundance than non-strange quarks. During hadronisation, numerous 
strange quarks can be bound with available non-strange quarks and strange 
hadron production will enhance. 

One can define a strangeness enhancement factor (E) as, 

1 dN"^^ 



1 dNPP \^-^) 

2 dy 

In Figj33l experimental data [108], [IDH] , in ^^^^=200 GeV Au+Au col- 
lisions for strangeness enhancement are shown. Data do show enhancement. 
meson (which is a hidden strange meson) production is also enhanced 
|110] . However, it is uncertain whether or not the enhancement is due to 
increased production in QGP or due to canonical suppression |112] |113] 
of strangeness in pp collisions. The canonical suppression arises from the 
need to conserve strangeness within a small, local volume, which limit the 
strangeness production in pp collisions relative to AA collisions. In the lan- 
guage of statistical mechanics, while canonical ensemble is applicable in pp 
collisions, in AA collisions, grand canonical ensemble is applicable. Assuming 
the systems correlation volume is proportional to Npart, canonical framework 
predict that yield per Npart increase with Npart as phase space restriction due 
to strangeness conservation is lifted. 
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9.3 Electromagnetic probes 

Photons and dileptons are considered to be important probes for QGP diag- 
nostic. They are called electromagnetic probe as they interact only electro- 
magnetically. Unlike the hadrons, which are emitted only from the freeze-out 
surface, photons and dileptons have large mean free path and are emitted 
from the entire volume. Total production is then obtained by convoluting 
their production rate over the 4- volume. 

9.3.1 Photons 

One of the problems with photon as a probe is the large background. For 
diagnostic purpose, one is interested only in 'direct photons', i.e. photons 
produced as a result of some collision process. However, QGP is a transient 
state and it ultimately transforms into hadrons. A large number of hadrons 
decay into photon. Decay photons constitute the back ground and needed 
to be eliminated to get the direct photon signal. Experimentally several 
methods have been devised to eliminate the decay photons, e.g. invariant 
mass analysis, mixed event analysis etc. However, due to overwhelmingly 
large number of background photons, the elimination could not be exact and 
experimental data on 'direct' photon production could not be obtained very 
accurately. There is definitely scope for further improvement. 

In a nucleus-nucleus collision, there are various sources of direct photons. 
They are briefly discussed below. 

(i) prompt photons: initial hard collisions produces prompt photons. 
Prompt photon production in a nucleon-nucleon reaction a + b — ?■ 7 + X 
can be calculated in leading order pQCD, 



= ^ y dxidxjfi/a{xi,Q'^)fj/b{xj,Q'^) 

x6{s + i + u)-^{ij ^-f + X) (9.10) 

TT dt 

where fi/a,b{x, Q^) is the parton distribution function, ^ is the elementary 
partonic cross section. Above equation scaled by the binary collision number 
gives the prompt photon production in nucleus-nucleus collision. 

(ii) Fragmentation photons: initial hard scattered partons may fragment 
into photon (g — )■ g + 7). 
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K ^ j dx^dxj—fi/a{xi,Q^)fj/bixj,Q'^)D^/k{z,Ql,) 

ijk=q,q,g 

x6{s + i + u)-i,{ij ^ k + l) (9.11) 
vr dt 

In Eq l9.11l -D-y/fc is the fragmentation function. As before, above equation 
should be scaled by the binary collision number to obtain the fragmentation 
photons in nucleus-nucleus collision. 

(iii) Pre-equilibrium photons: in nucleus- nucleus collisions, an extended 
QCD medium is produced. Before achieving local thermal equilibrium the 
medium is in pre-equilibrium stage. Photons will be emitted in the pre- 
equilibrium stage also. However, it is difficult to distinguish pre-equilibrium 
photons from thermal photons. Simulations with Parton Cascade Model |114] 
indicate that pre-equilibrium photon production equals the thermal photon 
at Pt=2 GeV. Low pj- photons are predominantly thermal. 

(iv) Thermal photons: photons, emitted from the (locally) equilibrated 
QGP and hadronic matter are called thermal photons. For QGP diagnos- 
tic purpose, these photons are most important. In the following, I will 
discuss briefly about thermal photons. For more detailed information, see 
[n3][TTB][TT7j . 

In the QGP phase, the important reactions for direct photons are; (i) 
annihilation: q + q ^ g + •y, (ii) Compton process: q + g ^ q + •j and (iii) 
bremsstrahlung: q + q ^ qq'j. 

In kinetic theory, photon production rate (per unit time per unit volume) 
from 1 + 2 — 3 + 7 process can be written as. 



E 



da 

d^p'^ 



Ui = Af [ ^P—^I^JJEi)fJE2)(27iy6\p'l+p^-p^.-p>') 

x|>(j2^^!^^!^[l±/3(^3)] (9.12) 

where Aii is the amplitude for one of the basic process, /(£')'s are the Fermi- 
Dirac or Bose distribution, as appropriate. The ± in the last term corre- 
sponds to Bose enhancement or Pauli blocking. Using the Mandelstam's 
variable, s,t,u, differential rate can be written. 
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X [1 ± fsiEi + E2- E)]e{Ei + E2- E)\laEl + hE^ + c 

(9.13) 



where, 



a = -{s + tf 

b = 2{s + t){Es- E2t) 

c = st{s + t)-{Es + E2tf (9.14) 

For massless particles, amplitude A4. is related to the differential cross 
section as, 

^-i^ (915) 

dt ~ 167rs2 ^^-^^^ 

For Compton and annihilation processes, the differential cross sections 

are, 



^^annhilation g^^^^ ^2 ^ ^2 



dt 9s^ ut 



dt us 



(9.17) 



A/" = 20 for annihilation process when summing over u and d quarks and 
A/" = ^ for the Compton process. 

Photon production rate from Compton and annihilation processes were 
first computed in |118]|119] |120] . Importance of Bremsstrahlung process was 
first considered in [121] |122] |123] , however, Landau-Pomeranchuk-Migdal (LPM) 
effect (when photon emission is suppressed due to multiple collisions) was 
neglected. Arnold, Moore and Yaffe |124] made a complete calculation in 
leading order. They also provided a simple parameterised form for easy use 
in hydrodynamics. Below, I list the results. 

Leading order photon emission rate from QGP: 
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= A{k)[\n{T/m^) + CUk/T)] (9.18) 



with, 



Ctot = \ H2k/T) + C2^2{k/T) + C,rem{k/T) + Canmh{k/T) (9.19) 

The leading log coefficient A{k) is given by, 

A{k) = 2a[dFY,q}P^nf{k) (9.20) 

K 

where nf{k) = [exp{k/T) + 1]^^ is the Fermi distribution function, and 
dp is the dimension of the quark representation {dp = Nc = 3 for QCD). 
gi=2/3 for up quark and -1/3 for down type quarks, moo = fi's^^/3 is thermal 
quark mass in the leading-order, Qs being the strong coupling constant, (a^ = 
gs/4:TT). For two flavor QCD, 

,,,, 407rT2 nf(k) 

A{k) = -^aas^f^ (9.21) 

C2^2{k/T), Ci,rem{k/T) and Cannih{k/T) all involve multidimensional in- 
tegrations, which can only be solved numerically. Numerical results for QCD 
plasmas are reproduced quite accurately by the approximate, phenomeno- 
logical fits [121] . 



C2^2{x) = 0.041x"^ -0.3615 + l.Ole"^-^^^', (9.22) 
[ iV7 0.548/o^(12.28+ 1/x) 

(9.23) 



y/l + x/lQ.27 

In the hadronic phase, photons are emitted in several reactions. Dominant 
channels are, (i)7r + 7r — p + 7, (ii) n + p — )■ 71 + 7. J. Kapusta, P. Lichard, and 
D. Seibert |118] |119] evaluated photon production rate from these channels. 
Role of Ai resonance in photon production was investigated in [125] . Later, 
in a more comprehensive analysis [126] , hadronic photon production rate in 
a meson gas consisting of light pseudo-scalar, vector and axial vector mesons 
{tt, K, p, K*, Ai) was obtained. A parameterised form was also provided. 
Below, I reproduce the parameterised reaction rates. 



104 




Figure 34: Invariant cross section (p + p) and invariant yield (Au + Au) of 
direct photons as a function of pt- See the text for details. 



Photon emission rate from Hadronic phase: 
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E' 



d^k 



F^{E)T'^-^ X exp 
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In the above equations, the E and T are in GeV and the rates are in unit 
of GeV~^fm~'^. The dipole form factor F{E) is, 



F{E) 



2A2 



2A2 



,A = IGeV 



(9.24) 



Photon production rate from the QGP phase and the hadronic phase has 
to be convoluted over the 4- volume to obtain total photon production during 
the evolution of the fireball, 



E- 



'd?k 



d'^x E 



,dR 

dFk 



J Tdrdxdydri (^E^^ (9.25) 



There are several simulations for direct photon production in relativistic 
heavy ion collisions. For demonstration purpose, I will show a few results. 
In FigjnH transverse momentum dependence of the invariant cross section 
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Fi gure 35: Thermal photon production in 0-10% Au+Au colhsions at a/s^^— 
200 GeV in different hydrodynamical models [130] [m] [l32] [133] [TBI] are 
compared with experimental data. 

in a/s^^=200 GeV p+p collisions and invariant yield in y^^jy=200 GeV 
Au+Au collisions are shown, p+p data are from |127] . The Au+Au data 
are from |128] |129] . The three curves on the p+p data represent NLO pQCD 
calculations. For pt > 2 GeV, the pQCD calculation is consistent with 
the p+p data within the theoretical uncertainties. The dashed curves in 
Au+Au data are obtained by scaling the photon yield in p+p collisions by 
the nuclear overlap function T^a- At low px <2.5 GeV, experimental Au + 
Au data are underpredicted. At low px, photon production increases faster 
than the binary NN collision scaled p + p cross section. The solid lines in 
Fig|33] are fits with an exponential plus binary collision scaled p+p data. 

In Figj35]experimental data on direct photon production in 0-10% Au+Au 
collisions are compared with different hydrodynamical model calculation. All 
the hydrodynamical simulations are comparable to the experimental data and 
with each others within a factor of 2. The results confirm the dominance of 
thermal radiation in the direct photon spectrum in low and intermediate pt 
range. 
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Figure 36: A schematic presentation of invariant mass dependence of dilepton 
production in high energy nuclear colhsions. 



9.3.2 Dileptons 

Similar to the photons, dileptons are also emitted at every stage of heavy 
ion collisions. In the QGP phase, a qq can interact to form a virtual photon, 
which subsequently decays in to a lepton pair or dileptons, g + g — > 7* — > /"'' + 
/~ (/=e or fi). In the hadronic phase, dileptons are produced from interactions 
of charged hadrons with their anti-particles e.g. 7r+ + vr" — p —)■/"'" + /~, 
from hadronic decays, e.g. 7r° —j-Z+H- /~ + 7, p —)■/"'" + /~, a; —)■ /"'"-I- /~, 
— >■ /"•" + Z~ etc. Dileptons are also produced in the Drell-Yan process (in 
the Drell-Yan process, a valence quark from the projectile hadron interacts 
with a sea anti-quark from target hadron to form a virtual photon, which 
then decays into a lepton pair). Depending on the dilepton invariant mass 
(M), dilepton emission can be classified into three distinct regimes. In Figl36l 
invariant mass (M) dependence of dilepton production in high energy nuclear 
collisions is shown schematically. One can distinguish three distinct regions, 
(i) M < M^{= 1.024 MeV) where dominating source of dilepton production 
is vector meson decays, (ii) < M < Mj/^ (= 3.1 GeV) where dominant 
source is the thermal production from QGP and (iii) M > Mj/^ dominated 
by primordial emission, decay of charmoniums etc. 

In the following, I briefly discuss dilepton emission rate in QGP and in 
hadronic resonance gas. As for photons, emission rates are to be convoluted 
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over the space-time volume to obtain production yield in nucleus-nucleus 
collisions. For more details see |115] . 

Dileptons production rate from a thermal system at temperature T, in a 
reaction, + a~ — )■ l'^ + l~ , where a is either a quark or pion, was obtained in 
|135]|136j . From QED, one calculate the cross section for e+-|-e~ — )■ fj,~^ + fi~, 



a{M) 





1 H 


2m2 




1 - 




3 M2 




V 

M2_ 






M2 



1/2 



(9.26) 



where M is the invariant mass of pair, a is the fine structure constant 

[a = 1/137) and mi is the mass of fi. For qq annihilation, the above is mul- 
tiplied by the color factor Nc=3, and factor reflecting the fractional charges 
of quarks. The modified cross section is. 



a(M) = Fga{M) 



iV,(2s + l)2 5^eJ, 
/ 



(9.27) 
(9.28) 



where s is the spin of quarks, e/ is the fractional charge and the sum is over 
the quarks flavors. In the hadronic phase, in the vector meson dominance 
mode, dilepton production goes via the reaction, 7r+ + tt" —)■ p —;■/+ + /~. 
The QED cross section is then multiplied by the Breit-Wigner form factor. 



150MeV 



[m 



Dilepton cross section from tt+tt annihilation then become, 

1/2 



a^{M) = F^{M)a{M) 



Ami 

TV 



(9.29) 



(9.30) 



Kinetic theory gives the reaction rate (number of reaction per unit time 
per unit volume) 



(27r)3 (27r)3 



/(Pl)/(P2) X Cr(a+a -)> /+/ )Vrel 



(9.31) 
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where, /(p) is the occupation probabihty at momentum p and energy E = 
Relative velocity Vrei can be computed as, 

Approximating /(p) = exp{—E/T), and integrating over five of the six 
variables, 

Ria+a- l^r) = — — - / a{z)z'^{z'^ - Azl)Ki{z)dz (9.33) 



(2vr) 

with z = M/T, Za = TTLa/T and Ki is the modified Bessel function of the 
first kind. 

Apart from the total number of lepton pairs emitted per unit space-time 
volume, R = j^, one is interested in several differential rates. They can 
be obtained from Eq J9.31l by appropriate change of variables. The rate for 
producing lepton pairs with invariant mass M is. 



^ 4m 



Production rate of leptons pairs with invariant mass M, momentum p and 
energy E (E = ^/p^TM^) can be written as, 



dN a(M) , . , 



^ Ami 



Presently, I will not discuss any hydrodynamic simulation for dilepton 
production. Dileptons have been measured in SPS energy {y/sj^j^=17.A 
GeV). Hydrodynamical simulations underpredict low invariant mass dilepton 
yield. In the low invariant mass region, dilepton emission is largely mediated 
by p meson. Now, properties of a hadron can change in a medium. Due to 
medium effect, p meson mass may drop, or its width increase. These two 
effects are considered as a likely cause of underestimation of low mass dilep- 
ton yield in hydrodynamic simulations. Currently, experimental data do not 
distingush between these two effects. For more detailed account, please see 
|115] . I may mention here that dileptons are possibly better suited to probe 
QGP than direct photons. In contrast to photons which are characterised 
by the transvesre momentum, dileptons are characterised by two parameters, 
the transverse momentum and invariant mass. For differential diagnosis, an 
increased degree of freedom may be useful. 
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Figure 37: (left panel) An experimental reconstruction of an event in pp 
collision. Two jet structure is clearly seen, (right panel) same for an event in 
Au+Au collisions. Notice that due to large number of particles in the final 
state jet structure is obscured. 

9.4 Jet quenching 

Jets are clusters of several hadrons, all moving in approximately the same 
direction. In hadron-hadron or in e~^e~ collisions one generally observe two 
jet structure with back to back momenta. In the left panel of FiglSTl an 
experimental reconstruction of final state particle trajectories is shown. The 
two jet structure is clearly evident in p+p collisions. In nucleus-nucleus 
collisions the jet structure is not obvious due to large number of particles 
(see the right panel of Fig J57|) . However, jets are there and one can use some 
'jet algorithm' to find them. 

Jet structure in hadron-hadron collisions can be understood qualitatively 
in perturbative QCD. The jet properties depend in general on two scales, the 
energy of the jet and its virtuality, specified by the largest possible transverse 
momentum of one of its subjets. The basic two body reaction 1 + 2 — )■ 3 + 4 
is in the partonic (quark or gluon) level. The scattered partons are highly 
virtual (squared 4-momentum transfer is large) and reduce their virtual- 
ity by radiating gluons or by splitting into quark-antiquark pairs. Such a 
parton branching is governed by Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
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Figure 38: Two particle correlation in the azimuthal plane, in A/i^^=200 
GeV, p+p, d+Au and Au+Au collisions is shown. Note the two peak struc- 
ture in P=P and d+Au collisions. In Au+Au collisions however, the away 
side peak is vanished. 



(DGLAP) equation. Finally the partons fragment into hadrons. The char- 
acteristic collimated hadrons from fragmentation of an outgoing parton are 
called jet. Naturally, the most common structure seen is the two jet event. 
Three jet events are also seen and results from reactions such as qq — )■ qqg. 

From theoretical consideration Xin-Nian Wang and Miklos Gyulassy pre- 
dicted the Jet quenching phenomena |137] . They argued that a partonic jet, 
if travel through a medium, will lose its energy by gluon emission. The energy 
degraded parton will ultimately fragment into less number of particles than 
it would have in absence of a medium. It can be demonstrated as follows: In 
leading order perturbative QCD, production cross-section for the hadron C 
in A + B^C + X reaction can be written as. 



xS{s + t + u) — ^{ab cd)Dc/c{zc,fJ'^) (9.36) 

In Eq J9.36l fa/A{x, Q"^) is the structure function of the parton a in the 
hadron A, fb/Bix, Q^) is the structure function of the parton b in the hadron 
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B. The is the factorisation scale. ^ is the partonic cross section for the 
reaction a + 6 — )■ c + rf. Dc/d^c, fJ''^) is the fragmentation function for the 
parton c into hadron C, fx"^ being the fragmentation scale, = ^ is the 
fraction of the partonic energy carried by the hadron C. K in Eg J9. 361 takes 
into account higher order effects. In nucleus-nucleus collisions, if colored 
medium is formed, the parton c can lose energy in the medium. If the parton 
c travels through a medium and lose energy A_Ec, then Zc = — )■ Zc{l — 
^^)~^. The fragmentation function Dc/c{z) is a rapidly falling function of 
z and increase in Zc will lead to reduced production for the hadron C. 

After the prediction of jet quenching phenomena, it was discovered in Rel- 
ativistic Heavy Ion Collider (RHIC) |138] |139] . In Figj38| di-hadron correla- 
tion in the azimuthal plane, in p-|-p, d+Au and Au+Au collisions are shown. 
The data are obtained in the following manner. A high pt trigger particle 
is fixed and in coincidence with the trigger particle, associated particles are 
measured as a function of the azimuthal angle. A0 = (^associate — 4>trigger 
is the difference of azimuthal angle between the trigger and associated par- 
ticles. In p-|-p and d+Au collisions, di-hadron correlations shows a double 
peak structure, which can be understood in terms of two jet events. The 
peak at A0=O is called the near side peak (nearer to the trigger particle) 
and the peak at A(j) = tt is called the away side peak (away from the trig- 
ger particle). In Au+Au collisions however, the away side peak is strongly 
suppressed. Strong suppression of the away side peak is the experimental ev- 
idence of jet quenching. The understanding is as follows: a di-jet is produced 
near the surface the medium. One of the jet escapes into the vacuum and 
fragments. The other enters the medium and loses its energy in the medium 
before fragmentation. 

Jet quenching lead to high pt suppression i.e. production of high pt parti- 
cles in A+A collision is less than that would have expected in a p-f-p collision, 
scaled by the collision number. High-p^ suppression is usually expressed in 
terms of the nuclear modification factor (Raa), 

dN^^/dy(PpT ,^ 
-^^^ ~ ^ — J pp/j M — ' (9.37) 

Taa in Eq J9.37l is the nuclear thickness function, calculable in a Glauber 
model. If AA collision is a superposition of pp collisions, the expected ratio 
is unity. Fig J39] shows the experimental measurements for Raa in a/s^^=200 
GeV Au+Au coUisions. Mesons are largely suppressed at high px- Photons 
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Figure 39: Nuclear modification factor in y/sj^j^=200 GeV Au+Au collisions 
for 7r°,7r^, t] and 7. The lines are theoretical calculations. 

are however are not suppressed. The observation is a strong evidence that 
high pt suppression is not an initial state effect, but a final state effect. High 
density colored medium is created in the collision and cause the suppression. 

Parton energy loss AE provides fundamental information on the medium 
properties. There are several models for the energy loss calculations, e.g. 
BDMPS [140] [141], GLV [142] [143]. I will not discuss them here. Just men- 
tion that energy loss in the medium is generally characterised by the following 
variables, 

(i) the mean free path A = 1/ (per), p the medium density and sigma the 
particle medium cross section, (ii) the opacity = L/X or the number of 
scattering centres in a medium of thickness L, (iii) the Debye mass mo{T) ~ 
gT, g the coupling parameter, (iv) the transport coefficient q = m'jy/X con- 
trols the radiative energy loss, responsible for jet quenching in the induced 
gluon bremsstrahlung picture, (v) the diffusion constant D = fiT character- 
ising the dynamics of heavy, non-relativistic particles, /i being the mobility 
of the particles, i.e the ratio of drift velocity and applied force. 

The various curves in Figj39] are from different theoretical calculations. 
The dashed curve shows a theoretical prediction using the GLV parton energy 
loss model |143j|144] . The model assumes an initial parton density dN/dy 
= 800 - 1100, which corresponds to an energy density of approximately 5-15 
GeV/fm3. The solid curves are predictions from reference [145] . The effect 
of parton energy loss was implemented through an effective modified frag- 
mentation function. The modified fragmentation function approximates the 
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medium effect in multiple parton scattering formalism. The dot-dashed curve 
in Figj39] is a theoretical result on Raa by considering only the collisional 
energy loss |146] . Theoretical predictions are approximately consistent with 
the experimental data. 

Our discussion on parton energy loss is rather sketchy. Interested reader 
may look into [147] for more detailed discussion. 

10 Summary 

In this short lecture course, I have discussed some aspects of relativistic heavy 
ion collisions. Our centre is actively engaged in experimental and theoretical 
studies on Quark-Gluon-Plasma. The emphasis of the course, naturally was 
directed to the study of QGP. I have discussed some topics at some length, 
some topics briefly and completely left out some topics. The choice of the 
topics is personal, which the author felt important enough for a student, 
pursuing his career in theoretical or experimental high energy nuclear physics 
should know. In future, I may extend the scope of the lecture note. I hope 
the students find this lecture course useful. I will be obliged to receive any 
comments or suggestions. Without any hesitation, reader may contact me 
by e-mail: akc@vecc.gov.in. 
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